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A  quadratic  cost  boundary  control  problem  for  a 

hyperbolic  system  on  a  domain  Q  ^         with  a  smooth  boundary 

r  is  discussed. 

Let  F  and  G  be  given  self-adjoint,  positive-definite 
bounded  operators  on  L^(q). 

Our  problem  is  to  minimize  the  functional  cost 
J(u,y)  =   |u|2q^^^^^  +  (i^'^y)C0,T]xQ       (y(T),Gy(T))^  over 
controls  u  €  L2([o,T]xr)  where  y  is  the  solution  to  the 
nonsymmetric  and  nondissipative  hyperbolic  system  of  first 
order : 

^  =  A(x,9  )y  in  [0,T]xi2 

'"ly  =  u  in  [0,T]xr 

y(0)  =  y^        in  Q 

with  y^  e  L^(Q) . 
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VJith  some  technical  assumptions  on  the  differential 
operator  A(x,d)   and  the  boundary  operator  M(x),  we  obtain  a 
semigroup  representation  for  the  hyperbolic  system. 

It  is  then  shown  that  our  problem  admits  a  unique 

optimal  control  u°  in  L^([0,T]xr)  with  the  corresponding 

trajectory  y^  in  C([0,T];   L'^(Q)).     VJe  show  that  the  optimal 

control  u°  can  be  expressed   in  a  feedback  form  via  a  Riccati 
2 

operator  on  L  (Q);  i.e., 

u°(t)   =  BP(t)y°{t)  a.e.   t  €  [0,T] 

where  B  is  an  unbounded  operator  from  L^(Q)    into  L^(r)  and 

P(t)   is  a  Riccati  operator  on  L^(Q)  which  satisfies  an 

appropriate  Riccati  equation. 

Under  some  smoothness  assumptions  on  F  and  G,  we  give 

additional  regularity  results  for  the  optimal  solution  and 

show  that  the  Riccati  operator  P(t)    is  the  unique  solution 

to  the  appropriate  Riccati  equations  within  a  certain  class 

of  one  paramater  families  of  operators  on  L^(Q). 

Finally,  it  is  shown  that  the  Riccati  operator  P(t)  for 
2 

the  L  -problem  in  which  F  and  G  are  assumed  to  have  no 
smoothing  properties  can  be  determined  as  the  strong  limit, 
uniformly  in  t,  of  a  sequence  of  solutions  to  the  Riccati 
equations  corresponding  to  regularized  problems  in  which  F 
and  G  are  assumed  to  have  some  smoothing  properties. 


CHAPTER  I 
INTRODUCTION 


1.1  Notations 


Let  Q   be  an  open  domain  in  R™  with  smooth  boundary. 

We  denote  by  c"(Q;r'^)   the  space  of  infinitely  dif f erentiable 

maps  from  Q   into  R*^  with  compact  supports,   and  by 

c"^j(Q,r'^)   the  restriction  of  C^(r"^;r'^)   to  Q. 
2  k 

Let  L  (Q;R  )   be  the  space  of  Lebesque  square  integrable  maps 

from  Q   into  R    with  natural  inner  product  (•,•)  and 

L  (Q.-R*^) 

norm     'I   ^         .    .     For  convenience  of  notation,  we  write 

2  11 
simply  L  (Q),    ('fO^  and   |.|^  if  there  is  no  confusion. 

Some  Sobolev  spaces  are  defined  as  follows:     For  s  >  0 

integer,   let  H^(Q;r'^)   be  the  completion  of  c"^^(Q,r'^)  with 

respect  to  the  norm 


lis  lis 

where  a  =   (a    ,...,a    )  multi-indices,    lal   =  a,   +...+  a  , 

-L  111  III  m 

a ,  a 

a  m 

a"'  =  -5   ^  

_    •  • .  . 

a  ^  a 

ox  ax 
1  m 

s  k  Ir 

Let  H^(Q;R  )  be  the  completion  of  c"(Q;r'^)  with  respect  to 
(1.1).     Similarly  Sobolev  spaces  of  noninteger  orders  are 


and 


1 


2 


defined  by  interpolation  as  in   [13].     As  a  matter  of 

convenience,  H^(Q;r'^)   and  H^(Q;r'^)     are  abbreviated  to 
s  s 

H  (Q)   and  H^(Q),     repsectively ,   if  no  confusion  arises. 

Let  J  be  an  interval  and  X  be  a  Hilbert  space.  We 
introduce  the  following  spaces: 


C(J;X)    is  the  space  of  continuous  maps  from  J  into  X,  with 
supremum  norm, 

L  (J;X)   is  the  space  of  measurable  and  essentially  bounded 

maps  from  J   into  X,   with  essential  supremum  norm, 

L^(J;X)    is  the  space  of  Lebesque  integrable  maps  from  J  into 

X,  with  norm   |y|  =  /    |y(t)|  dt, 

L-^(J;X)  ^  ^ 

2 

L  (J;X)    is  the  space  of  Lebesque  square   integrable  maps  from 

J  into  X,   with  norm   lyl^  =  /    |y(t)|^dt,  and 

L^( J;X)  ^ 

h1(J;X)    =   {y|y  €   l'^{J;X)    and  ^  y  €   L^iJjX)}    with  norm  i 
H   (J;X)  L^(J;X)        '^^       L^{J',X)  ■ 


In  view  of  Fubini's  theorem,  we  have  L^(J;L^(Q))  = 

2 

L  (JxQ).     For  a  bounded  linear  operator  P  from  a  Banach 

space  X  into  a  Banach  space  Y,  we  denote  bv  IIPII 

^  X-*-Y 

I    f"P     I^^Iy  operator  norm.     If  there  is  no  confusion, 

r  Ix"^ 

we  write  simply  ii  ph  .  vje  also  denote  by  X'  the  dual  space  of 
a  Banach  space  X. 
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1.2     Problem  Formualtion  and  Preliminary  Material 

Let  Q   be  an  open  bounded  domain  in  r'^  witn  smooth 
boundary  r . 

VVe  consider  a  first  order  differential  operator  of  the  form 

m 

(1.2)        A(x,9  )    =     Z     A  .  (x)d  .    +  C(x)  , 

j  =  l      J  J 

where  9^  =  5/5 x^  and  the  coefficients  A^ ( x )   and  C(x)  are 

smooth  n  X  n  matrix-valued  functions  defined  on  Q. 

We  make  the  following  assumptions  on  the  operator  A(x,5): 

(A.l)      (Strict  Hyperbolicity )  i 

For  all  nonzero  vectors  (u  )    in  r"'  and  all  x  in  Q,  the 

m  1  m 

matrices     s     A.(x)w-  have  n  distinct  real  eigenvalues. 

M 

For  X  in  r,   let  A  (x)   =     E     A.(x)n.(x)  where 

j  =  l     ^  ^ 

(n^(x) , . . . ,n^(x) )   is  the  inward  unit  mormal  to  the  boundary 
r  at  x. 

(A. 2)     (Noncharacteristic  Boundary) 

The  matrix  A^{x)   is  nonsingular  for  each  x  in  r. 

Remark.     From  the  continuity  of  the  A j  '  s  and 
assumptions   (A.l),   (A. 2),  we  deduce  that  A^(x)  has  n 
distinct  real  eigenvalues  different  from  zero  for  each 
X  in  r.     Henceforth,  we  assume,  without  loss  of  generality, 
that  A^(x)   has  k  distinct  negative  eigenvalues  and  (n-k) 
distinct  positive  eigenvalues  for  each  x. 
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VJe  impose  a  boundary  condition  on  the  operator  A(x,d) 
with  the  aid  of  a  boundary  operator  M(x)   which  is  a  smooth  k 
X  n  matrix-valued  function  defined  on  r  of  rank  k  and  annihi- 
lates all  vectors  in  the  null  space  of  a'^(x)    for  each  x  €  r. 

We  now  consider  the  following  initial  boundary  value 
problem:   For  T  >  0,   f  e  L^([0,T]    x  Q),   u  €  L^([0,T]    x  D  and 
y^  €  L^(Q ) , 

(1.3)       a  ^y  =  A(x,d)y  +  f  in    [0,T]    x  Q 

My  =  u  in   [0,T]    x  r 

y(0)   =  y^  in  Q 

Let  us  denote  Q  =   [0,T]   x  Q  and  E  =   [0,T]    x  r.     ■  ' 

We  recall  the  definition  of  strong  solutions  to  (1.3). 

Definition.     A  function  y  in  L^(0)    is  said  to  be  a 
strong  solution  to  (1.3)    if  there  exists  a  sequence  {y  }  in 
C^qj(Q)   and  v  in  L  (Z)   such  that  Mv  =  u,  ^  . 

I^n  -  "  " 

bn  -  He  0 

htYn  -  A(x,5)y^  "   ^Iq  ^  0 
ai^d  l^n^^^   "  ^oIq       0  as  n  ». 

It  is  well  known   [15]    that  by  using  a  partition  of  • 
unity  and  change  of  local  coordinates,   the  above  mixed  prob- 
lem (1.3)   can  be  transformed  into  a  problem  on  a  half-space. 
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Thus  we  assume,  without  loss  of  generality,   that  Q  is 
the  half-space  {x|x  =   ( x^^ ,      ,  .  .  .  ,  x^)   €  r"^,   x^  >   O}  with 
r  =  {x|x  €  r"',   x^  =  0}  .     For  the  half-space  Q,   the  matrix 
a'^(x)    in  (A. 2)   becomes  Aj^(x)    for  each  x  in  r.     In  view  of 
finite  speed  of  propagation  in  hyperbolic  systems,  we 
require,  without  loss  of  generality,   that  the  matrices  A j ' s , 
C  and  M  are  constant  outside  a  compact  set  in  Q  [15]. 

From  the  assumptions   (A.l)   and   (A. 2),   by  a  smooth 
change  of  coordinates,  A-^  can  be  written  in  the  following 
standard  form: 


(1.4) 


A^(x)  = 


A^  (x) 


0 


a|(x) 


for  X  e  Q  where 


A^(x)  = 


with  negative  funtions  a^,...,aj^, 


and  A^ ( X )  = 


'k+1 


'k+1 


with  positive  functions 


For  simplicity,   we  use  the  following  notation: 
X  =   ^^^1:^2' *  •  • ^   (Xj^,x''")  and 

y  =  (yi'-'-Vkiyk+i'-- wYn)  =  (y~'  y^). 

As  to  the  boundary  operator  M  for  the  half-space  Q,  we 
simplify  as  follows.     Let  us  partition  M  as   [HiN]    where  H  is 
a  k  X  k  matrix  and  N  is  a  k  x   (n-k)  matrix. 


with  A-[^(x)    in  the  standard  form  (1.4),   it  is  well  known 
that  a  necessary  condition  for  (1.3)   to  be  well-posed  is 
that  H  be  nonsingular   [5,   15]  .     Since  M  and  H'-'-M  impose 
essentially  the  same  boundary  conditions,  we  assume  that 

• 

(1.5)       M  =   [IlN]   where  I  is  the  k  x  k  identity  matrix. 

In  order  to  make  the  mixed  problem  (1.3)  well-posed  in 
the  Kreiss'   sense,  we  need  to  introduce  the  so  called 
"Kreiss'   condition"   formulated   in  terms  of  eigenvalues  and 
generalized  eigenvalues  as  follows. 

For  a  moment,   we  assume  that  all  the  coefficient 
matrices  A j ' s  and  M  in  (1.3)   are  constant,  and  we  then 
consider  the  boundary  value  problem: 


We  take  the  Fourier  and  Laplace  transforms  of  (1.6) 
with  respect  to  the  tangential  variables  x"*-  =  (x2,...,x^) 
and  the  time  variable  t  with  the  dual  variables  cj  = 
(u2  <■  •  •  •  f Uj^)     and  s  =  r  +  iq,  respectively.     We  then  obtain 


m 


(1.6)       ^  ^.Y  = 


S  A -9  .y  in  0 
j  =  l     J  ^ 


My  =  0 


in  Z 


A 


(1.7) 


dy 


A^      (sl  - 


i     E     A  .(jj  .  )y 
j  =  2     J  J 


for  x 


>  0 


1 


(1.8) 


A 

My 


0 


for  X,   =  0 


A 


where  y  is  the  Fourier  and  Laplace  transforms  of  y. 
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We  recall  the  following  definition  from   [8] . 

Definition .     A  complex  number  s  =  r  +  iq  is  said  to  be 
an  eigenvalue  of  (1.3)    if  there  exists  a  nontrivial  solution 
<t)(x^)   of  the  ordinary  differential  equation  (1.7)   and  (1.8) 
in  L^(  [0,=»)  ;r")  . 

-1  . 

Since  A,    (si  -  i     E     A.oo.)   has  exactly  k  eigenvalues 

j  =  2         ■'  V 

with  negative  real  parts  for  Re  s  >  0   [8] ,   it  is  well 

known  that  any  solution  ^  of  (1.7)   in  L^([0,«))   can  be 

expressed  as  a  linear  combination  of  k  linearly  independent 

normalized  solutions  (j>  ^  (  s  ;.),...,$  ^  (  s  ;•  )   in  L^([0,oo)). 

Hence,   from  (1.8),  a  number  s  with  positive  real  part  is 

k 

an  eigenvalue  of   (1.3)    if     z     m  0).(s;.)C,  =  0  for  some 

j=l         J  J  ,  • 

nonzero  vector  (c^,...,c^);   i.e.,   the  square  matrix 

(M())^  ,  (s;  •),...  ,M(|)j^(  s; .)  )      is  singular.  -  .  . 

Definition.     A  pure  imaginary  number  iq  is  said  to  be  a 
generalized  eigenvalue  of  (1.3)    if  the  matrix 

(M(l)^(s;.  )  ,  .  .  .  ,M(l>j^(s;.  )  )   approaches  a  singular  matrix  as 
r  >  O"*",  where  s  =  r  +  iq. 

We  assume  the  following  condition. 
(A. 3)      (Kreiss'  condition) 

For  each  boundary  point   (0,x^)    in  r,   the  constant  coefficient 
problem  (1.3)   that  arises  by  freezing  the  coefficients  A's 
and  M  at  (0,xM   has  no  eigenvalues  or  generalized 
eigenvalues  with  nonnegative  real  parts. 
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We  refer  to   [8,   15,   18]    for  other  equivalent  conditions. 

Now  we  are  in  a  position  to  recall  basic  results  on 
well-posedness  and  regularity  for  (1.3)   from  [15]. 

Theorem  1.1     VJell-posedness . 

Assume   (A.1)~(A.3).  •  • 

Then,   for  f  €       {Q)  ,         €  L^{Q)   and  u  €  h'^iZ),   the  variable 
coefficient  mixed  problem  (1.3)   admits  a  unique  strong 
solution  y  in  L^(Q).     Furthermore,   the  following  estimate 
holds : 


(1.9)        |y(t)|     +  |y 


+ 

0 


Y 


4  Iq  +    Ms  +    IYqI^]    for  0  <   t  <  T, 

where         is  a  constant   independent  on  f,  y^,   u  and  t. 

Theorem  1.2  Regularity. 
Assume   (A.l )  ~   (A.3) . 

Let  s  be  any  positive  integer.     Suppose  that  f  e  H^(0), 

Yq  €  H^(Q)   and  u  €  H^(i:)   satisfy  the  following  compatibility 

conditions: 

^t  ^lt=0  "  °  "lt=0  "  °  integers  0  <   j   <  s . 

Then  the  solution  y  to  (1.3)   has  the  following  regularity: 

y  e  H=(0),   y(t)   e  H^(S2)   and  y|^  e  H^il),   0  <   t  <  T. 
Moreover  for  0  <   t  <  T, 


(1.10)       y(t)  1    o  +  IyI    .^  +  |y| 

s     I    ls,Q       Mols,Q       I  \s,Z^ 
where  Cg  is  a  constant  independent  on  f,  y^,   u  and  t. 

With  theorem  1.1  and  1.2  at  hand,  we  are  now  in  a 
position  to  formulate  our  optimal  control  problem. 
We  introduce  the  state  equation; 

(1-11)  a       =  A{x,d )y  in  0  . 

My  =  u  in  E 

y(0)  =  Yq  in  Q 

where  A(x,5)   and  M  are  assumed  as  in  theorem  1.1, 
y^  €  L^(Q  )   and  u  €  L^(E  )  . 

From  theorem  1.1,   for  any  function   (control )   u  in 

2 

L  (Z),     there  corresponds  a  unique  solution   ( trajectory)  to 
(1.11)   which  will  be  denoted  by  y'^.  - 

F  and  G  be  any  bounded,  self-adjoint  and  positive- 
definite  operators  on  L^(Q). 

We  consider  the  following  functional  cost: 

(1.12)       J(u,y^)   =  1    [|u|2  ^   (y^Fy")^  -H   ( y^  T )  ,Gy^  T )  )  1 


Q 


2 

for  u  in  L  ( E )  . 


We  note,   from  theorem  1.1,   that   (1.12)    is  well-defined. 
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Remark .     A  more  general  functional  cost  such  as 

J(u,y")  =  \  +  (y"-y^F(y^V))Q  +  ( y^  ( T) -y^  ( T)  ,  ' 

G(y'^(T)-y^(T)  )  )^]  ,  where  y^  is  a  given  function  in  I?  {Q)  with 

y  (T)  €  L  (Q),  can  be  treated  in  the  study  which  follows. 
However,  for  the  sake  of  simplicity  of  notation,  we  shall 
restrict  ourselves  in  what  follows  to  functional  cost  (1.12). 


Optimal  Control  Problem; 

(CP.)     Minimize  the  functional  cost  J(u,y^)   over  u  in 

2 

L  (2)     subject  to  the  state  equation  (1.11). 
Our  basic  questions  are 

(i)  Does  the  (CP.)   admit  a  unique  optimal  control 
in  L^(E)?  < 

(ii)  If  yes  in  (i),   can  the  optimal  control  be 
expressed  in  a  Riccati  feedback  form? 

Our  aim  is  to  give  affirmative  answers  to  the  above 
questions.     To  elaborate  (ii),  our  goal  in  this  work  is  to 
express  the  optimal  control  u°  in  a  feedback  form; 

u°(t)   =  BP(t)y°(t),  a.e.   t  €  [0,T] 

where  B  is  an  unbounded  operator  from  L^(Q)   into  L^(r),  and 
P(t)   is  a  Riccati  operator  on       (Q)  which  satisfies  a 
corresponding  Riccati  equation. 

More  details  are  presented  in  sections  1.3  and  1.4. 


11 


1.3     Survey  of  Literature 


A  numDer  of  authors  have  studied  optimal  control 
problems  for  dynamics  which  include   first  order  hyperbolic 
systems.     Here  we  briefly  mention  some  of  the  results 
already  obtained  in  this  field  and  compare  them  with  the 
results  obtained  in  this  thesis. 

In  the  fundamental  treatise  by  Lions    [12]    on  optimal 
control  problems  for  partial  differential  equations, 
boundary  control  problems  for  second  order  hyperbolic 
equations   (not  for  first  order  systems)   are  briefly 
mentioned  under  the  assumption  of  smooth  controls. 

Russell   [16]   has  studied  boundary  control  problems  for 
first  order  symmetric  hyperbolic  systems  in  one  spatial 
variable,   through  the  method  of  characteristics.     It  is  well 
known  that  the  method  of  characteristics  cannot  be  applied 
to  problems  in  several  spatial  variables.        v.-;  - 

Curtain  and  Pritchard   [2,   3]    have  studied  optimal 
control  problems  describee  by  evolution  equations.  When 
they  specialize  their  results  to  boundary  control  problems 
for  first  order  hyperbolic  systems,   it  is  needed  to  assume 

that  the  control  spaces  are  l2([0,T];   H^Cr))    (rather  than 
L^( [0,T] ;    L2(r) ) ) . 

The  work  that  is  closest  in  spirit  to  ours  is  a  paper 
by  Vinter  and  Johnson   [18],    in  which  they  discuss  Riccati 
equations  for  hyperbolic  systems  with  boundary  controls. 
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A  brief  summary  of  the  results   in   [18]    is  provided  here. 
Vinter  and  Johnson   [18]    consider  functional  cost   (1.12)  , 
under  the  assumption  F  =  0,   that  is, 

(1.13)  J(u,y'')   =    |u|2  +   (y^(T),  Gy^d))^. 

By  using  a  variational  techniques  as   in   [12] ,  they 
characterize  the  optimal  control  in  terms  of  the  solution  z 
of  the  two-point  boundary  value  problem; 

(1.14)  d  ^.z  =  -A*z  -  Fy  in  0 

M*z  =0  in  Z  •  ' 

z(T)   =  Gy(T)  in  9  ^ 


and 


a  j.y  =  A(x,a  )y  in  0 

My  =  u  in  2 

Y(0)   =  Yq  in  Q 


where  A*  is  the  formal  adjoint  of  the  operator  A(x,d)  and  M 
is  the  adjoint  boundary  operator  of  M. 

Orily  under  the  supplemental  assumption  that  G 

"'^P^  ^^^^^  into  H^(Q),   they  derive  an  appropriate  Riccati 

differential  equation,   and  then  synthesize  the  optimal 
control  in  a  feedback  form.     The  main  technical  tools  in 
[18],   besides  variational  techniques,   are  the  regularity     "  ■ 
theorems  from  [15],   together  with  trace  theorems. 

To  summarize,   the  following  comments  pertaining  to  [18] 
are  in  order: 
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(i)  Operator  G  in   (1.13)    is  assumed  to  carry  L  (Q) 
into  H^(Q).     This  means  that  the  penalization  of 
the  terminal  state  y(T)   occurs   in  a  space  of 
distributions,   instead  of  L^-space. 

(ii)  Functional  cost   (1.13)   does  not  penalize  the 
trajectory,   i.e.  does  not  contain  the  term 
(y,Fy)Q  of  (1.12). 

However,   the  penalization  of  the  trajectory  is 
important,  because  it  is  precisely  responsible  for  the 
minimization  of  the  whole  energy  of  the  system,   which  in 
turn  is  a  primal  task  of  engineering  implementations.  The 
techniques  of   [18]   cannot  be  extended  to  treat  the  important 
case  when  the  trajectory  or  the  terminal  state  is  penalized 
in  the  physically  meaningful  L^-spaces.     This  is  due  to  the 
fact  that  the  regularity  theorems  in   [15]   are  not  applicable 
to  the  case  of  F      o.     To  understand  this  point  better,  let 
us  mention  the  argument  of   [18]   which  relies  on  differen- 
tiability of  the  solution  z  of  (1.14). 

In  order  to  deduce  dif f erentiabilty  of  the  solution  z 
by  applying  Rauch's  regularity  theorems,  some  compatibility 
conditions  are  needed  for  the  terms  Fy  and  Gy(T)   in  (1.14). 
But  without  assuming  F  =  0,   the  compatibilty  conditions 
cannot  be  satisfied.     For  this  reason,  one  should  use 
different  techniques  from   [18],    in  order  to  overcome  the 
limitations  discussed  in  (i)   and  (ii). 

Motivated  by  the  physical  significance  of  the  problem, 
we  study  Riccati  feedback  synthesis  with  boundary  controls 
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in  L  -spaces;    i.e.,   the  trajectory,   the  terminal  state  and 
the  control  are  all  penalized  in  L^- topolog ies .     This   is  a 
distinguishing   feature  of  our  work   from  the  existing 
literature.   In  this  sense,  we  extend  and  complement  [18]. 

More  precisely,  we  do  obtain  feedback  realization  of 
the  optimal  control  as  postulated  in  (CP.),  without  any 
smoothness  assumptions  on  the  operators  F  and  G.     It   is  also 
shown  that  the  feedback  operator  P(t)  does  satisfy  a  Riccati 
integral  equation.     Our  techniques  different  from  those  in 
[18]    are  based  on  more  refined  analysis  of  regularity  of  the 
solution  to  the  mixed  problem  (1.3).     This  analysis  leads  to 
a  more  accurate  description  of  the  trace  of  the  Riccati 
operator  P(t),  which  is  a  necessary  ingredient  in  the 
derivation  of  Riccati  equations. 

In  this  analysis,  we  use  a  combination  of  partial 
differential  equation-  and  semigroup-  methods  together  with 
so-called  "direct  approach"  as  in   [9,10].     It  should  be 
pointed  out  that  penalizations  in  L^-spaces,   instead  of 
distributions,   lead  to  major  mathematical  difficulties 
related  to  the  trace  of  the  Riccati  operator  P(t).  - 

The  Riccati  differential  equation  for  the  original 
L  -problem  may  not  have  meaning,  even  in  the  distributional 
sense.     Here  we  use  the  regularization  technique  to  overcome 
the  difficulties. 

In  a  regularized  problem,  we  show  that  the  Riccati 
operator  P(t)   is  the  unique  solution  to  appropriate  Riccati 
equations  in  a  certain  class  of  operators  on  L^(Q).  The 
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original  L  -problem  is  then  recovered  as  a  limit  of  the 
regularized  problems.     This  way,  we  show  that  the  Riccati 
operator  for  the  original  problem  is  determined  as  a  strong 
limit  of  Riccati  operators  for  a  sequence  of  regularized 
problems. 

From  the  numerical  point  of  view,   the  regular i zat ion 
results  can  be  considered  as  a  tool  to  construct  a 
convergent  algorithm  for  computing  the  Riccati  operator  of 
the  original  L^-problem.     Finally,   the  question  of 
regularity  for  the  optimal  solutions   to  our  control  problem 
is  answered  under  some  assumption  on  F  and  G.  Precise 
statements  of  our  results  are  given  in  the  next  section. 


1.4     Main  Results 

In  order  to  state  our  main  results,  we  need  some 
preliminary  material. 

Let  A  be  the  operator  defined  by 

(1.15)       Ay  =  A(x,d)y  for  y  in  the  domain  D(A) 

where  D(A)   =  {y  €  L^{Q)  \   Ay   in  L^{Q)   and  My  |     =  0}  . 
It  is  well  known  [15,   18]    that  A  generates  a  strongly 
continuous  semigroup  S(t)   on  L^(Q). 

VJe  next  introduce  the  "Dirichlet"  map  "D"   (a  natural 
extension  from  the  boundary  r  into  the  interior  Q)  defined 
by  "Du  =  y"  where 
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(1.16)  A( x,d )y  =  ky  in  Q 
My  I p  =  u  i  n  r 

for  some  sufficiently  large  constant  k  >  0  (Dependence  of  D 
on  k  will  not  be  explicitly  indicated).     We  show  in  section 

2.1  that  the  above  Dirichlet  map  D  is  a  bounded  linear 

2  2 
operator  from  L  (D    into  L  (Q)    (theorem  2.1). 

With  the  help  of  the  operators  A  and  D,   we  give  an 

abstract  semigroup  model  for  the  solution  of  (1.9). 

In  fact,  we  define  an  operator  L  as 

(1.17)  (Lu)(t)   =   (A-kl)         S( t-s)Du(s)ds  for  0  <   t  <  T 

and  u  in  the  domain  of  L  which  is  a  dense  subspace  of  L^(Z). 
Then  we  prove,   in  section  2.2,   the  following  theorem. 

Theorem  2.2     Semigroup  representation. 

(i)  The  operator  L  is  a  linear  bounded  operator  from 
L^(Z )    into  C(  [0,T]  ;   L^(Q)  )  . 

(ii)  The  solution  y  to  the  mixed  problem  (1.9)    is  given 
by  the  semigroup  representation:  '- 

y(t)   =  S(t)y^  -   (Lu)(t)    for  0  <   t  <  T. 

In  section  3.1,  we  show  that  our  control  problem  (CP.) 
admits  a  unique  optimal  control  u°  in  L^(Z). 

Now  we  are  in  a  position  to  state  our  main  results 
(Theorem  3.2,   4.1,   5.1,   5.2,   and  5.3). 
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Theorem  3.2     Regularity  of  optimal  solution. 

Under  the  assumption  that  F  and  G  are  compact  operators 

s  1 
on  H  (Q)   for  all  0  <  s  <  ^-z   it  is  shown  that 

( i)  u°  is  in  H^(Z ) 

( ii)  y°  is  in  H^(0) 

for  any  initial  data         in  H^(Q). 

Let  us  denote  by  $ ( t ,  s )  ,   0  <   s  <   t  <  T,   the  evolution 
2 

operator  on  L  (Q)   which  describes  the  evolution  of  the 
optimal  trajectory  y°. 

Theorem  4.1     Feedback  synthesis  of  optimal  control. 

The  optimal  control  u°  can  be  expressed  as  follows. 

(i)  u°(t)   =  BP{t)y°(t)    for  t  in  [0,T], 

where  B  =  D*(A*-kI)   which   is  given  by  By  =  A~y~ | ^ 
for  y   in  D( A* ) . 

(ii)  P(t)    is  a  bounded  operator  on  L^(Q),  satisfying 
the  following  Riccati  integral  equation  (the  first 
Riccati  equation): 

For  x,y  in  L^(Q),    t   in  [0,T], 

(I  R.E.)        (P(t)x,y)Q  =  /J  ($(s,t)x,F$(s,t)y)^ds 

+  (BP(s)$(s,t)x,BP(s)$(s,t)y)^ds 
+  ($(T,t)x,G$(T,t)y)^. 

As  a  corollary  of   (I  R.E.),   we  obtain  that  P(t)    is  self- 
adjoint  and  positive-definite. 
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Control  problems  with  smooth  observations  are  treated 
in  theorems  5.1  and  5.2. 


Theorem  5.1     Riccati  equations. 

Under  the  assumption  that  F  and  G  map  L^(Q)    into  the 
domain  D(A*),  P(t)   satisfies  the  following  Riccati  equations 

(i)  (Riccati   integral  equation,   which  we  call  the 

second  Riccati  equation): 
2 

For  x,yinL(Q),   0<  t<T, 

(II  R.E.)        (P(t)x,y)^   =  /J   (S(s-t)x,FS(s-t)Y)^ds  - 

l\  (BP(s)S(s-t)x,BP(s)S(s-t)Y)^ds  + 
(S(T-t)x,GS(T-t)y)^  . 

(ii)  (Riccati  differential  equation): 
For  x,y  in  D(A)   and  a.e.   t  in  [0,T], 

(D.R.E.)       ^   (P(t)x,y)^   =  -^""'^y^Q   -  (P(t)Ax,y)^ 
-   (P(t)x,Ay)^  +  (BP(t)x,BP(t)y)^ 
with  terminal  condition  P(T)  =  G. 

Theorem  5.2  Uniqueness. 

Under  the  same  assumptions  as   in  theorem  5.1, 
P(t)    is  the  unique  solution  to  the  Riccati  equations  in 
theorem  5.1  within  the  class  P  ,  where  P     is  the  class  of 
one  parameter  families  {Ht)}^^^  of  self-adjoint,  positive- 
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definite  and  bounded  operators  on  L  (Q)   such  that 

6P(»)    is  bounded   from  L^(Q)    into  l'°([0,T];  L'^(r)) 
and     BP(«)S(.)B*   is  bounded  from  L^(r)    into  L^(E). 

Let  P^  be  the  correspond i ny  Riccati  operators   for  the 
control  problems  with         and  G^,  in  the  functional  cost  J(u, 
Finally,   we  obtain  the  following  convergence. 

Theorem  5.3  Convergence. 

Let  F  and  G  be  any  self-adjoint,   positive-definite  and 
bounded  operators  on  L  (Q). 

Suppose  that  F^  and  G^  map  L^(Q)    into  D(A*),   and  that 

>  F  and        >  G  strongly  in  L^(Q)   as  n  >  ». 
Then,   Pp,(t)  -»•  P(t)   strongly  in  L'^iQ)   as  n  ->  =  uniformly  in 
Moreover,   for  x,y  in  L^(Q)   and  uniformly  in  t, 

(P(t)x,y)^   =  j'^   (S(s-t)x,FS(s-t)y)^ds  +   (  S  (  T- 1 )  x  ,GS  (  T- 1 )  y  )  ^ 
-  lim         (BP   {s)S(s-t)x,BP   (s)S(s-t)y)  ds. 

Remark .     The  significance  of  theorem  5.3   is  two  fold: 

(i)  it  serves  as  a  tool  for  proving  theorem  4.1  (ii), 

(ii)  it  enables  to  prove  convergence  of  an  effective 
numerical  algorithm  to  determine  the  feedback 
operator  P(t)    for  the  L^-problem. 

The  rest  of  the  paper  is  devoted  to  proofs  of  the  main 
results  as  follows. 
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In  section  2.1,  we  prove  the  existence  and  regularity 
of  the   "Dircihlet  map"  D  (theorem  2.1).     Section  2.2  gives 
the  proof  of  theorem  2.2.     In  section  3.1,  we  explicitly 
characterize  the  optimal  control  u°  and  optimal  trajectory 
y°  (theorem  3.1).     In  section  3.2,   we  give  regularity 
results  for  the  optimal  solutions   (theorem  3.2  and  theorem 
3.3).     In  section  4.1,   the  evolution  operator  is 
introduced   (lemma  4.2).     Section  4.2  gives  the  proof  of 
theorem  4.1   (i)   and  a  partial  proof  of  second  part   (ii).  In 
section  5.1  and  5.2,   we  give  the  proof  of   theorem  5.1  and 
5.2.     Finally,    in  section  5.3,  we  give  the  proof  of  theorem 
5.3  and  we  complete  the  proof  of   theorem  4.1   (ii).  Also, 
examples  of  sequences  {f^}   and  {g^^}    in  theorem  5.3  are 
given. 


CHAPTER  II 

SEMIGROUP   REPRESENTATION   FOR   HYPERBOLIC  SYSTEMS 


2.1     Dirichlet  Map 

In  order  to  introduce  a  semigroup  representation  for 
the  solution  to  (1.11),   we  need  to  extend  functions  defined 
on  r   to  functions  on  Q . 

An  operator  which  maps  boundary  functions  into 
specially  defined  interior  extensions  is  said  to  be  a 
"Dirichlet  map." 

In  this  section,  we  show  that  the  Dirichlet  map  "D" 


defined  in   (1.16)    is  a  bounded  linear  operator  from  L^(r 


2 

into  L  (Q) . 


We  consider  the  following  boundary  value  problem: 

(1.16)       A(x,d )y  =  ky  in  Q 
My  =  u  in  r 

for  u  in  L2(r)  where  k  is  constant  >  0.     The  existence  of 
our  Dirichlet  map  D  follows  from  the  next  theorem. 

Theorem  2.1. 

There  exists  a  constant  k  >  0  large  enough,   such  that 
for  any  u  €  L^U),    (1.16)   has  a  unique  solution  y  in  L^(^). 
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Furthermore,   the  following  energy  inequality  holds: 

(2.1)  UlQ+|y|p<c|u|^ 

for  a  constant  c  independent  on  u. 
The  proof  will  be  given  later. 

Corollary.     The  Dirichlet  rrap  D  is  a  bounded  linear 
operator  from  L^(r)    into  L^(Q). 

Remark.     The  above  result  was  first  proven  in   [11]  for 
symmetric  problems  with  dissipative  boundary  conditions,  and 
then  was  extended  to  symmetrizable  problems  with  dissipative 
boundary  conditions  in   [4].     m  our  nonsymmetric  and 
nondissipative  case,   the  proof  of  theorem  2.1  will  be  given 
by  using   "Kreiss'   symmmetrizer  technique"  as  in  [8]. 

Proof  of  Theorem  vve  may  assume,  without   loss  of 

generality,   that  C(x)   =  0   identically.     We  shall  first  prove 
the  energy  inequality   (2.1)   and  then  the  existence  of  the 
solution  to  (1.16).     suppose  for  a  moment  that  the  coeffi- 
cient matrices  A  -  s  and  M  are  constant  frozen  at  a  boundary 
point   (0,xl)    in  r.     we  apply  the  Fourier  transform  to  (1.16) 
for  the  tangential  variables  x^  with  the  dual  variables 
=   (0)2,  .  .  .  ,a)j^)  .     Then   (1.16)  becomes 


(0 

A 


^^•^^       h  =    tkl-i     z     A  (o.]y  for  x,    >  0 

1  j=2      J    J  1 

A  A 

My  =  u  c 

for  x^  =  0. 
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Since  A]^   is  nonsingular,    (2.2)  becomes 
dy  -1  "'a 

^  =  A       [kl-i     E     A.oj.Jy  for  x,    >  0 

A  A 

My  =  u  for  x^  =  0 

which  is  the  same  sort  of  problem  Kreiss  studied  in  [8]. 

Now  we  may  rewrite   (2.2)    in  pseudo-differential   form  as 

'       ^^'^^       ^^^^  ^  =  B(x,a))i   ,  > 


0 


1    A  A 

'M(x    )y  =   u,  =  0 


where  B(x,a,)  is  a  symbol  of  the  pseudo-differential  operator 
(of  first  order)   corresponding   to  the  differential  operator 


m 

kl  -     E     A  .  (x)d  .  . 
j  =  2     J  J 


Since  A-^   is  nonsingular,   we  have 

''^'^^       ^  =  V(x,a))y   ,   x^   >  0 

A  A 

My  =  u,  =  0 

where  V(x,co)  =  A~^  (  x )  B  (  x  ,a) )  . 

Now  we  use  the  Kreiss  symmetrizer  which  is  constructed 
locally  in  conical  neighborhoods  of   (x,(o,0)  [8]. 

VJe  recall  the  following  results  from   [8].  , 

Lemma  2.1.  There  exists  a  symmetrizer  for  (2.4)  which 
is  a  pseudo-differential  operator  (of  order  zero)  whose  , 
symbol  R(x,a;,k)     with  L2-norm  independent  of  k  satisfies 
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A 

( i )       R  is  Hermitian, 


A 


(ii)     R  is  homogeneous  of  degree  zero  in   (k,oj)  for 
2       I    I  2 

k     +   |w|     >   1  and  is  a  smooth  function  of  all  its 

variables  and  of  the  coefficient  matrices  A^'s  ana 

^  I    1  2         I    ,  2 

(ill)    (y,Ry)   >  5|y|     -  e|u|      for  all  vector  y  such  tnat 

My  =  u  in  r . 

A 

(iv)     Re  RV  >   6k  where  5   >  0  and  s   >  0  are  constants 
independent  on  k,aj  and  x. 

Now  we  are  in  a  position  to  prove  the  energy  inequality 

A 

(2.1).     Multiplying  the  symbol  R  to   (2.4)   and  taking  the 

A 

inner  product  with  y,  we  obtain 


A 

^    d\7  ^  AAA 


'1 

It  is  understood  that  the  inner  products  in  (2.5)  are  taken 
for  the  variables   { ,^ ^ , . , . ,^   ) .     From  (iv),    (2.5)  becomes 


(2.6)       Re   (R^,   J)^^   >  6k|;|2^ 

On  the  other  hand,   using   integration  by  parts  with 
respect  to        and  the  fact  that  =  0,  we  obtain 

A 


A 


A       A  J 


From  ( i ) ,  we  have 

A 


(2.7)       Re   (R  liL_,   y)     =iRe[-(tRt)  U^ 
dx^'  -^'q       2     ^   ^Y'^y^p  -  (y,  y)g] 
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A  A 


dR 

(2.7)  becomes 


From  (iii)  and  the  fact  that   (y,  2£_  ^3        ^^^^  ^raer, 


for  some  constant  c  >  0   independent  on  y  ana  u. 
By   (2.6)   and   (2.8),   we  obtain 

(2.9)        (6k-c)    |;|^^  ^  1  5|J|2   <   1  ,|i|2. 
Since  the  Fourier  transform  is  unitary,   we  may  rewrite   (2.9)  as 

i 

By  choosing  k  large  enough  so  that  6k-c  >  0,   we  have 

^   for  some  constant  c. 

From  (2.11),   the  energy  inequality   (2.1)   can  be  easily 
derived . 

Next  we  turn  to  the  proof  of  existence. 
First  we  notice  that  since  the  operator  A  generates  a 
strongly  continuous  semigroup  on  L^{Q),    there  exists  a  real 
number  r  such  that  the  set  {^|Re  X  >  r}    is  contained  in  the 
resolvent  set  of  A.     Hereafter,  we  choose  and  fix  k  large 
enough  to  be  sure  that   (2.11)   holds  and  k  >  r. 
Then   (A-kl)   has  a  bounded  inverse  on  jJiQ).     That  is, 
(2.12)        [A(x,& )-kI]y  =   f   in  Q 
My  =  0  in  r 

has  a  unique  solution  in  l2(q)    for  any  f   in  L^{Q). 


|y|p  <  c 

1  ^  1 
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To  continue  with  our  proof,  we  first  consider  smooth 
boundary  functions  u  in  (1.16),    i.e.   u  €  H^(r).     It   is  well 
known  that  there  exists  a  smooth  extension  u  in  H^(Q)   of  u. 

Let  u  =   (u   ,u"^)   where  u~  =  u  and  u"^  =  (0,...,U). 
Then  u  e  H^(Q;r")   and  A(x,3)u  €  L^(Q). 
Thus,    from  (2.12),    it  is  easily  seen  that 

(2.13)        [A(x,5)-kl]y  =  -A(x,5)u  +  ku   in  Q 
My  =  0  in  r 

admits  a  unique  solution  y  in  L^(Q). 

From  (1.5),  Mu  =  u  in  r.     Let  y  =  y  +  u.    ;  ■  ' 

Then  from  (2.13),  y  is  a  solution  to 

(1.16)       A(x,a  )y  =  ky   in  S2 
My  =  u  in  r. 

This  implies  the  existence  of  solution  to  (1.16)  for 
smooth  boundary  data  u.     By  using  a  standard  density  argument, 
we  extend  our  results  to  L^-boundary  data  u  as  follows. 
For  a  given  u  in       (T ) ,   we  choose  a  sequence  {u^}^   in  H^(r) 
such  that  u^  .  u  in  L^(r)   as  n  .  ».     gy  energy  inequality 
(2.1),   the  solution  to  (1.16)    is  unique   if   it  exists.  Let 
yn  be  the  solution  to  (1.16)   for  the  boundary  data  u„;  i.e., 

A(x,5 )y^  -  ky^  =  0   in  Q 
Myn  =  %  in  r. 

Then  by   (2.1),   for  any  positive   integers  m,  n. 
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(2.14)  y  -y    L   +    I V  -V    I      <   c    lu  -u  I 

Since  {u^}^  is  a  Cauchy  sequence  in  L^(r),    from  (2.14) 

we  deduce  that  {y^}^  and   {y^lpl^^  are  Cauchy  sequences  in 
2  2 

L  (Q)  and  L  (D  respectively.  Thus,  there  exist  the  limits 
y  in  L^(Q)   and  y^   in  L'^(r)   such  that 

(2.15)  kn'^la  ^  °  l^n"^rlr  ^  0  as  n  >  =0. 
From  (2.15),  we  obtain  My^  =  u. 

This  means  that  y   is  a  strong  solution  to  (1.16). 
Proof  of  theorem  2.1   is  thus  completed.  , 


Remarl<_.     if  the  coefficients  Aj(x)   in  A(x,d)  are 
symmetric  and  block  diagonal,   then  we  can  obtain  (2.1)  with 


the  aid  of  a  simple  symmetrizer  R  = 


el^  0 


0 


"n-k 


where  I, 


and  are  k  X  k  and   (n-k)   x   (n-k)    identity  matrices 

respectively,   and  e   is  a  small  positive  number. 


That  is,  A.  = 
3 


^  _ 

A  . 

: 

0 


0 


-  A^  where  A~  are  k  x  k  symmetric 


matrix  functions  and  A^  are  (n-k)  x  (n-k)  symmetric  matrix 
functions  for  j=l, . . . ,m. 

Suppose  that  y  is  a  solution  to  (1.16). 

Using  integration  by  parts,    it  is  easily  seen  that 


(2.16)        (RAi5^y,y)^   =  -   1    (RA^y,y)^   -   ^  (Ry,(5^A^)y) 


Q 


and  (RA.9^y,y)^ 


^  (Ry,(a jA^)y)^   for  j   >  2. 
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Since   (RA(x,5)y  -  R1<y,y)^  =  0,  we  have,   from  (2.16),  that 


ra  _ 


(2.17)  (RA^y,y)^  +  ( Ry , (d ^ ) y ) ^  +  2k(Ry,y) 

Boundary  condition  My  =  u  on  r  and   (2.17)  imply 

(2.18)  e(A^u,u)^  -  2e(A^u,Ny^)^  +  e  ( A~Ny"^  ;  Ny^  )  ^ 


Q 


+  +     +  ™ 


+   (A^Y    .y    )p   +   2k(Ry,y)      +     z      (Ry,(d.A^)y)      =  0. 

J  =  ]^  J    J  " 


Using  the  fact  that  A^   <  0 ,  a|  >   0,   and  that  A.,   N  are 
constant  for   |x|   large,  we  obtain  that  with  e  >  0  small  and 
k  >  0  large. 


|y|p  <  c 

1  ^  1 

This  completes  the  proof  of   (2.1)    for  this  special  case. 

In  the  sequel,  we  need  a  representation  of  the  trace 
operator  in  terms  of  the  adjoint  of  the  Dirichlet  map  D. 

This  can  be  done  as  follows. 
Let  A*  be  the  formal  adjoint  of  A(x,5). 
By  simple  computation,  we  have 


m 

^^•2°^                =  -.^^   A^(x)a.y  -  a.(A^(x)y)    for  y  €  D(A^ 

where  the  domain  D(A*)   =  {y   |   A*y  g  l2(q),  M*y|^  =  o},  and 

the  adjoint  boundary  operator  M*   is  given  by 
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(2.21)        M*  =    [-(At)    ^   n'^A,    :    I     ,  ] 

1  1    .  n-k 


Lemma  2.2.     The  operator  D*(A*-kI)    is  given  by 

D*(A*-kI)y  =  A^  y|~  for  y  €  D(A*). 

Rroof_.     For  y  €  D(A'')   and  u  €  L^(r),   by  Green's 
formula,  we  have  v 

(2.22)  (A*y,Du)^   =   (  y ,  A(  x  ,d  )  Du )  ^  +  (A^y,Du)j,. 

Thus,   by   (1.16),    for  y  G  D(A*)    and  u  €  L^(r),    (2.22)  becomes 

(2.23)  (A*y,Du)^   =   (y,kDu)^  +  (A^y,Du)^. 
Since  y  €  D(A*),   M*y  =  0  on  r .     From   (2.21),   we  have 

(2.24)  -(a|)~^  N^A^y"  +  y^  =  0   on  r . 

On  the  other  hand,   from  (1.16),   MDu  =  u  on  r,  i.e., 
(Du)~  +  N(Du)"^  =  u  on  r . 

Thus,   from  (2.23)   and   (2.24),  we  obtain 
'      (2.25)       (  (A*-kI)y,Du)2  "   ^^\y~'^)^  for  all  u  in  L^(r). 

This  completes  the  proof  of  the  lemma. 

Remark.     From  Lemma  2.2,    it   is  clear  that   the  operator 
D*(A*-kI)    is  densely  defined;    i.e.,   D(A*)    is  contained  in 
the  domain  of  DMA*-kI).     Thus,   it  has  the  adjoint. 
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2.2     Input-Solution  Formula  for  Hyperbolic  Systems 

In  this  section,  we  give  a  semigroup  representation  for 
the  solution  to  (1.11);  i.e.,  we  give  the  proof  of  theorem  2 
We  recall  the  operator  L  defined  in  section  1.4, 

(1.17)        (Lu)(t)   =   (A-kl)         S( t-s)Du(s)ds,   U   <   t  <  T 

for  u  in  the  domain  of  L. 

We  first  observe  that  L  is  closed.     This  can  be  seen 
[7]    by  viewing  L  as  a  product  of  the  operator   (A-kl)  which 
has  a  bounded  inverse,  with  the  closed  operator 
S( t-s)Du( s)ds  which  is  actually  bounded.  , 

We  next  show  that  L  is  densely  defined  on  L^(S). 
To  do  this,  we  compute  A        S(t-s)f(s)ds  for  f  €  H^([0,T]; 
D(A)).     Then  we  obtain,   after  integrating  by  parts, 

(2.26)  A        S(t-s)f(s)ds  = 

S^'^-s)         f(s)ds  +  S(t)f  (0)-f  (t)  . 

The  right  hand  side  of   (2.26)   can  be  extended  to  all  f  in 
H^( [0,T] ;   L^iQ)) .  Thus, 

(2.27)  A  Jl  S(t-s)Du(s)ds  =  /J  s(t-s)D  ^  u ( s ) ds 

+  S( t)Du(0)-Du( t) 

is  well  defined  for  u  6  H^([0,T];  L^(r)). 
Consequently,   L  is  densely  defined  on  L^iz). 
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In  order  to  prove  theorem  2.2,  we  need  the  following 
preliminary  result. 

Lemma  2.3.     The  operator  D*A*S*(»)    is  bounded 
from  L^(S)    into  L^(E).     That  is, 


!l   |D*A*S*(t)x|2  dt  <   c    |x|2   for  x  €  L^(Q) 
where  c  is  a  constant. 

£roof_.     Me  consider  the  following  mixed  adjoint  problem 
backward  in  time:   for  x  e  L^(Q)   and  f  €  L^(Q), 

(2.28)  &       =  -A*z-f   in  0 
M*z  =  0  in  Z 

Z(T)    =  X  in  Q  - 

It  is  well-known   [15]    that  the  mixed  problem  (1.3) 
satisfies  the  Kreiss'   condition  if  and  only  if  its  adjoint 
problem  satisfies  the  Kreiss'   condition  with  time  traveling 
backward.     Thus,   Rauch's  results  can  be  applied  to  (2.28). 
By  theorem  1.1,   the  following  estimate  holds: 

(2.29)  |z|q  +   |z|^  <   ^M^lg  ^   I'^IqI    ^'^^^  ^  constant  c. 

On  the  other  hand,   the  solution  z  can  be  expressed  in 
terms  of  the  adjoint  semigroup  S*(t)  as 

(2.30)  z(t)   =  S*(T-t)x  +  Jl  S*(s-t)f(s)ds,   0  <   t  <  T. 
Since,   by  theorem  1.1,   z  has  a  well-defined  trace  on  r  and 
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M*z|j,  =  0,   lemma  2.2  can  be  applied  to  (2.30).     That  is, 
D*(A*-kI)z( t)   =  A~(Z(t)|^)~  for  a.e.   t.  Thus, 

(2.31)  |D*(A*-kI)z|     =    |a"(zL  )"L   <   c   IzL  L 

for  a  constant  c  >  0,  since  A~  is  constant  outside  a  compact 
set.     From  (2.29),    (2.30)   and   (2.31),   we  obtain 

(2.32)  |D*(A*-kI)S*(T-. )x  + 

D*(A*-kI)   /y  -S*(s-. )f (s)ds|      <   c   [|fL   +   Ixl  ] 

for  f  e  L^(0)   and  X  e  L'^(Q).  „       .  j. 

By  taking  f  =  0   in   (2.32),   we  have    '        ■  -i*-- 

(2.33)  |D*(A*-kI)sMT-.  )x|^  <   c   |x|^   for  x  €  L2(Q). 

Since  kD*S*(T-.)    is  bounded  from  L^(Q)   into  L^(Z),  from 
(2.33),   so  is  D*A*S*(T-.).     Thus,   lemma  2.3   is  proved. 

Now  we  continue  with  the  proof  of  theorem  2.2. 

Proof  of  Theorem  2.2   ( i) .     By  taking  x  =  0   in  (2.32), 
we  obtain 

(2.34)  |D*(A*-kI)  /'[.)S*(s-.)f(s)ds|j, 

<   c   I f 1^  for  f  €  L^(0) . 
For  a  moment,   let  us  define  an  operator  L  as 
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(2.35)        (Lf)(t)   =  D*(A*-kI)    JI  S*(s-t)f (s)ds 
for  f  e  l'^(O)  . 

Then  by   (2.34),   L  is  a  bounded  linear  operator  from  L^iQ) 

2  — 
into  L  (Z).     It  is  easily  seen  that  the  adjoint  of  L  as  an 

2  2 

operator  from  L  (Q)   into  L  (E)   coincides  with  the  operator  L 
defined  in  (1.17).     Since  L  is  closed  and  densely  defined, 
L*  =  L;    i.e.,   L  is  the  adjoint  of  L. 
Thus,   L  is  bounded  from  L^(S)    into  L^(0). 

Next  we  prove  that  L  is  bounded   from  L^(E)  into 
l"( [0,T] ;    L^(Q ) ) . 

To  do  this,    for  u  e  L^(E)    and   f  €  L^([0,T];   L^(Q)),  we 
compute   (Lu,f)g.     By  lemma  2.3  and  Schwarz's  inequality, 

(2.36)    (Lu,f)^  =  ((Lu)(t),f(t))^^dt 

"  /n  /o   (u(s)  ,D*(A*-kI)SMt-s)f  (t)  )„dsdt 


1 

^  /o   f/ol"(^Hr^s]^    [/o|D*(A*-kI)S*(t-s)f(t)  l^ds] 


1 

^o   l^'lz    t/oh*^^*-'^^)S*(s)f  (t)  |2ds]2dt 


<    |u|^  |D*(A*-kI)S*(.)f(t)  Ldt 

^   l^lr    1^1   1  9  3  constant  c  >  0 

lU  [0,T]  ;   L^(Q)  ) 

That   is,    for  u  6  L^(S)   and   f  6  L^([0,T];  L^{Q)), 

au,f)g  <   c   |u|      |f|  .     j„  f^^^^   Lu  is   in  the 

L   ( [0,T] ;    L^(Q) ) 
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dual  space  of  L^([0,T];   L^(Q)),   which  in  turn  implies  that 

Lu  is   in  L°°([0,T];   L'^(Q))    for  each  u  6  l.'^iZ). 

In  order  to  prove   that  Lu   is   in  C([0,T];   L2(Q))  for 
2 

u  €  L   (E),   we  use  a  density  argument. 

From  (2.27),    for  u  €  H^([0,T];   L2(r)),     we  have 

(2.37)        (Lu)(t)    =         S(t-s)D  ^  u(s)ds  +  S(t)Du(0) 

-  Du(t)   -  k  /q  S( t-s)Du( s)ds. 

It  is  easily  seen  that  Lu   is   in  C([0,T];   l''^(Q))  for 
u  €  H^(  [0,T]  ;    L^(r)  )  .  ..  •• 

Since  H^([0,T];   L^(r))    is  dense  in  L^(E),   for  any  u  in 
2 

L  (Z),   there  exists  a  sequence   {u^}^  in 
H-'-([0,T];   L^(r))    such  that    lu  -u  I  0   as   n  ->•  =. 

O  2., 

For  any  positive  integers  m,n. 


^"^"'^"n  o  =   Ilu  -Lu  I 

"  C([0,T];    L^(Q))        '  "  l"([0,T];  L^(Q)) 


<  c     u  -u 
I  m 


^Ij,   for  some  constant  c  >  0. 


Hence  {Lu^}^  is  a  Cauchy  sequence   in  C([0,T];  L^(Q)). 

By  completeness  of  C([0,T];  L^(Q)),   there  exists  the  limit 

in  C([0,T];   L^(Q))    such  that 

|Lu  -y|  0  as  n  ->■  CO. 

c( [0,T] ;    L^(Q) ) 

On  the  other  hand,   {Lu^}^  converges  to  Lu  in  l'"([0,T]; 
L  (Q))   since  L  is  bounded  from  L   (E)    into  l"([0,T];  L^(Q)). 
Thus,   Lu  =  y  which  completes  the  proof  of  part  (i). 
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Proof  o£  theorem  2.2   (ii).     V>/e  first  consider  a  stnooth 
control  u  €  H-"-  (  E  )  . 

Suppose  that  y  is  the  solution  to  (1.11)  with  u  e  H^CE). 

Let  y  =  y  -  Du.     Then,   it   is  clear  that  y  is  the  solution  to 

(2.38)       9 ^y  =  A(x,5)y  "  D  ^  u  +  kDu  in  Q 

My  =0  in  Z 

y(0)   =  y^  -  Du(0)  in  Q. 


That  is, 

(2.39)  5^y=Ay-D^u+  kDu  in  Q 

y(0)   =  y^  -  Du(0) .  in  Q 

Using  the  variation  of  parameter  formula,  we  obtain,  from 
(2.39),  that 

(2.40)  y(t)   =  S( t) [yQ-du(O) ] 

-         S(t-s)[D  ^  u(s)-kDu(s)]ds 

=  S(t)y^-S(t)Du(U)   -         S(t-s)    D  ^  u(s)ds 

+  k  S(t-s)Du(s)ds. 

Since  y(t)  =  y{t)  +  Du ( t ) ,  from  (2.27)  and  (2.40),  we  obtain, 
for  u  e  H-'-  (Z  )  , 

(2.41)  y(t)   =  S(t)y^  -   (A-kl)  S ( t-s ) Du ( s ) ds . 
That   is,   y(t)   =  S(t)y^  -   (Lu)(t)    for  u  €  H^(Z). 
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Since  H   (E)    is  dense   in  L   (E)   and  L  is  boundea   in  L   ( >] )  ,  the 
formula  (2.41)   can  be  extended  to  all  of  L^(E)   by  density- 
argument. 

This  completes  the  proof  of  theorem  2.2. 
From  theorem  2.2   (i),  we  have  the  following. 

Corollary  2.2.     For  each  t  €   [0,T],  let 
l'^u  =   (Lu)  (  t)  ,   u  €  L^(Z  )  . 

Then  L*-  is  a  bounded  linear  operator  from       ( E )    into  L^(Q) 
uniformly  in  t,  0  <   t  <  T. 


CHAPTER  III 

CHARACTERIZATION  OF  OPTIMAL  CONTROL  AND  REGULARITY 
3.1     Characterization  of  Optimal  Control 

In  this  section,  we  first  establish  the  existence  and 
uniqueness  of  the  optimal  control  for  our  control  problem 
with  the  aid  of  theorem  2.2.  .:  ^ 

We  then  provide  an  explicit  representation  of  the 
optimal  control  and  the  corresponding  trajectory  in  terms  of 
the  intial  data.     This  will  be  done  by  using  the  Lagrange 
multiplier  method  [14]. 

In  order  to  prove  the  existence  and  uniqueness  of  the 
optimal  control,  we  recall  the  following  result  from  [12]. 

Lemma  3.1.     Let  VI  be  a  real  Hilbert  space.  Suppose 
that  Tt(»,«)   is  a  continuous  symmetric  and  coercive  bilinear 
form  on  W. 

Then  for  any  continuous  linear  functional  ^  on  v/,  the 
functional  cost  J(u)   =  7i(u,u)   +  ^{u)   admits  a  unique  element 
u°  in  a  closed  convex  subset        of  W  such  that  Jiu^)   <  J(u) 

for  any  u  €  W  ,   u  it  u  . 

o  o 

From  lemma  3.1,  we  obtain  the  following  lemma. 
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Lemma  3.2.     Our  control  problem  (CP)   admits  a  unique 
optimal  control  u°  in  L^(Z). 

Proof_.     By  using  the  semigroup  representation  in 
theorem  2.2,  we  may  rewrite  J(u,y)  as 

(3.1)        J(u,y)    =  i   |u|^  +  ^   (S( . )yQ-Lu,FS( . )yQ-FLu)^ 

+  I  (S(T)y^-L^u,GS(T)y^-GL^u)2 
=  ■rt(u,u)   +  (j)(u)   +  c  for  u  €  L'^(E) 
where  ii(u,v)   =  i  (u,v)^  +  i  (Lu,FLv)^^  +  i  (l'^u,GL^v)^  , 
<t>(u)   =   (Lu,FS(.  )yQ)Q  +   (L^u,GS(T)y^)^  and 
c  =  I  (S(.)y^,FS(.)y^)Q  +  i  (S(T)y^,GS(T)y^)^. 

Then  it  is  easily  seen,   from  theorem  2.2    (i)   and  Corollary 
2.2,   that  71  and  ^   satisfy  all  the  requirements  of  lemma  3.1. 
This  proves  the  lemma. 

Our  next  task   is  to  give  an  explicit  representation  of 
the  optimal  solution.     This  is  done  in  the  next  theorem. 

Theorem  3.1. 

The  optimal  control  u°  and  corresponding  trajectory  y° 
are  given  by  the  following  formulae:     for  y^  €  L^(Q), 

(i)  u°  =    (I  +  L*FL+L^*GL^)"^    [L*FS(-)    +  L'^*GS(T)]y^ 
in  L^(Z  )  , 

(ii)  y°  =  S(.)y^-Lu°  in  C([0,T];  L^(Q)). 
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Proo_^.     By  theorem  2.2   (ii),   we  may  reformulate  our 
control  problem  as  follows: 

(3.2)     Minimize  the  functional  cost 

j(u,y)  =  1  |u|2  +  1  (y,Fy)^  +  i  ( S ( T ) y^-L^u ,GS ( T ) y^-GL^u ) ^ 
subject  to  the  constraints 


(3.3)       y  -  S(»)y^  +  Lu  =  0  over  u  €  L^(E) 


o 


by 


At  this  point,   we  introduce  Lagrange  functional  defined 

(3.4)         i:(u,y,z)   =  J(u,y)   +   (  z  ,  y-s  (  •  )  y^+Lu  )  ^ 

for  u  e  L^iZ),  y  e  L'^(O)   and  z  €  L'^(O). 

It  is  clear  that  the  Lagrange  functional    C  is 
continuously  Frechet  dif ferentiable  with  respect  to  u,  y  and 
z   in  the  corresponding  spaces  L^(E)   and  L^(0).     Since  the 
optimal  control  exists  uniquely   in  L^(E),    it  can  be 
determined  as  a  critical  point  of  the  Lagrangean    C    [14] , 
That  is,   if  let     C  ^  and  be  the  Frechet  derivatives 

of     C    with  repsect  to  u  and  y  respectively,   then  the 
optimality  conditions  are  given  by 

(3.5)  £    (u°,y°,z°)   =  0  and     C   (u°,y°,z°)   =  0. 

y 

Thus,  by  simple  computations,    (3.5)    is  reduced  to 

(3.6)  z°  +  Fy°  =  0 

u°  +  L*z°  -  L^*GS(T)y^  +   L^*Gl'^u°  =  0. 
From  (3.6),  we  obtain 
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(3.7)  u°  =  L*Fy°  +  l'^*G(S(T)y^-L^u°). 
Since  y°  =  S{»)y^  -  Lu°,    (3.7)  becomes 

(3.8)  (I  +  L*FL+l'^*GL^)u°  =   L*FS(.)Yq  +  l'^*GS(T)Yq. 

We  notice  in   (3.8)    that  since  F  and  G  are  self-adjoint  and 
positive-definite,   so  is  the  operator  L*FL+l'^*Gl''- ,   as  an 
operator  on  L  (Z).     Thus,   the  spectrum  of  L*FL+l'^*Gl'^  is 
contained  in  the  nonnegative  real  axis  [17]. 
Therefore,  -1   is   in  the  resolvent  set  of  L* FL+l'^*GL"'"  . 
This  implies  that  I+L*FL+L^*GL^  has  a  bounded  inverse  on 
L  (Z).     Since  the  right  hand  side  of   (3.8)    is   in  L   (Z),  we 
obtain  the  first  part  of  the  theorem. 

The  second  part   (ii)   follows  from  the  first  part  (i) 
and  theorem  2.2  (ii). 

Thus,   the  theorem  is  proved. 

3.2     Regularity  of  Optimal  Controls 

The  purpose  of  this  section  is  to  studY  the  regularity 
of  optimal  control   (theorem  3.2).     In  the  general  case, 
without  assuming  smoothness  on  the  operators  F  and  G,   we  can 
only  say  that  u°  is   in  L  (Z).     If  we,   however,    impose  very 
mild  assumptions   (compactness  type)   on  F  and  G,   it  can  be  said 
that  the  optimal  control  u°  is,    in  fact,   smoother  than  . 


41 


This  result,   besides  being   important   in   its  own  right, 
has  also  important  application  to  numerical  analysis  of  the 
control  problem  where  regularity  of  optimal  control  plays  a 
crucial  role  in  determining  rate  of  convergence  of 
approximation. 

In  order  to  prove  theorem  3.2,   we  assume  additionally 

that 

(A. 4)     F  and  G  are  compact  operators  on  H^(Q)   for  all 
0  <  s  <  -J. 

We  define   the  following  spaces;    for  nonnegative 

integers  s,     '  ' 
o 

(3.9)  H^{Z) 
T 

H^(Q) 

Then  we  have  the  following  lemma. 

Lemma  3.3.     VJith  s  =  0  or  1 ,    the  following  operators 
are  bounded: 

o 

(i)  L:   H^(S)  H^(0) 

T     °s  s 

(ii)  L   :   H   (S)   >  H  (Q) 

(iii)  S(t):   hJ(Q)       H^(Q)   uniformly  in  t 

( iv)  S( • ) :   H^(Q )   >  H^(0) 

T 

(v)  L*:   H^(0)   -V  R^(Z) 

(vi)  l'^*:   H^(a  )   ■>   H^(i:  )  . 


=  {u  I  u  e  H^(Z ) ,  u I =  0}  and 
=  {y   I   y  e  H^(0)  ,  y|^^^  =  0} . 
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2£Oof_.     We  recall  the  mixed  problem  (1.11)   and  an 
adjoint  problem: 

5  ^.y  =  A(  x,d  )  Y  in  0 

iMy  =  u  in  S 

y(0)  =  yQ         in  Q 

5        =  -A*z-f     in  Q 

M*z  =  0  in  Z 

Z{t)   =  in  Q  f 

L^(Q),   f  e  L^(0)   and  u  €  L^(S).     Then,   as  in 
Chapter  II,   the  solution  y  and  z  are  given  by  the  following 
semigroup  formulae; 

(3.11)  y(t)   =  S(t)y^  -   (Lu) (t) 

(3.12)  z(t)   =  S*(T-t)z^  +         S*(s-t)f (s)ds. 

We  now  apply  the  regularity  theorem  1.2  to  (3.11)  and 

(3.12).     Then  by  taking  y^  =  0  in  (3.11),  we  obtain  (i)  and 

(ii).     Similarly,  by  taking  u  =  0  in  (3.11),  (iii)  and  (iv) 

are  obtained. 

Let  Zq  =  0   in   (3.12).     From  theorem  1.2, 
T 

z(t)  =        S*(s-t)f (s)ds  has  a  well-defined  trace  on  the 

T 

boundary  in  h^(E)   if  f  e  H^(0).     Thus,  we  obtai  n  ( V )  by 
using  lemma  2.2  and  theorem  1.2  together  with  the  fact  that 
A^(x)    is  smooth  on  r  and  constant  for   Ixl  large. 


and 


(3.10) 


for 


y  ,  z  € 

^  o  o 
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Similarly,   by  takiny  f  =  0   in   (3.12),    (vi)   can  be 
obtained. 

Proof  of  theorem  3.2.     By  interpolating  between  s  =  U 
and  1   [13],   lemma  3.3  holds  for  all  0  <  s  <  1,   s  it  -j. 
On  the  other  hand,   for  smooth  domain  Q,  the  following 
identifications  hold   [13]:   for  0  <   s  <  -j, 

o 

H    (S)  =  H^(Z) 
T 

H^(0)  =  H^(0)  .  - 

H^(Q)  =  H^(Q) .  ' 

By  using   (3.13),  we  may  rewrite  lemma  3.3   as  follows. 
For  0  <   s  <        the  following  operators  are  bounded: 

(  i)  '  L:   H^(2)   >  H^(0) 

(ii)  '  l'^:   ti^(l)  ->■  H^(Q) 

(iii)  '  S(t):   H''^(Q)  >  H^(Q)   uniformly  in  t 
(  iv)  •  S(»  )  :  H^(Q  )  H^(0) 

(V)  •  L*:   H^(0)  H^(i:) 

(vi)  •  l"^*:   h^(Q  )       H^(E  )  . 

Hence,   for  0  <   s  <  ^,  L*FL+l'^*Gl'^  is  a  bounded  operator 

on  H   (E).     Since  F  and  G  are  compact,   so  is  the  operator 
T  T 

L*FL+L  *GL   ,   being  a  sum  of  compositions  of  compact 

operators  with  bounded  operators.     Thus,    in  order  to  show 
T  T 

that  I+L*FL+L  *GL     has  a  bounded  inverse  on  H^(S),   by  the 


(3.13) 
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Fredholm  alternative,   it  is  sufficient  to  prove  that  -1  is 

not  an  eigenvalue  of  L* FL+L ■'^*Gl'^  on  H^iZ).     On  the  other 

hand,   we  know  in  the  previous  section  that  -1   is  not  an 

eigenvalue  of  L*FL+L^*Gl'^  on       ( E )  .     Since  H^(£)   is  contained 

in  L^(E),     -1  cannot  be  an  eigenvalue  of  L*FL+l'^*GL^  on  H^{Z}. 

Our  next  step  is  to  show  that  L*FS(»)y  +L'''*GS(T)y  is 

o  o 

in  H^(Z)   for        €  H^(Q).     This  follows  immediately  from 
( iii)  '   -   ( vi) •  . 

Thus,   u°  =   (  I+L*FL+l'^*Gl'^)~-'- [L*FS(  .  )y^+L'^*GS(T)y^]  is 
in  H^(E )    for  y^  €  H^(Q )  . 

Part  (ii)   follows  from  part   (i)   and  theorem  3.1  (ii). 
Thus,   the  theorem  is  proved. 

Remark .     If  s  >         then  the  identifications  (3.13)  do 

not  hold.     Thus,    in  view  of  theorem  1.2,   the  operator 
T  T 

L*FL+L  *GL     is,   in  general,   not  bounded  from  H^(E)   into  itself 

In  the  sequel,  we  need  also  regularity  of  optimal 
control  for  regularized  problem.     More  precisely,  we  have 
the  following  result. 


Theorem  3.3. 

Assume  that   (A. 5)    F  and  G  map  L^(Q)    into  hJ^(Q). 
Then  u°  e  H^(S)    for  any  initial  data  y^  e  L^(Q). 

^roof_.     By   (A. 5),  ql'^u  and  GS(T)y^  are  in  h-'-(Q)  any 
u  6  L^(S )   and  y^  €  L^(Q )  . 
From  lemma  3.3   (vi),   we  deduce  that 
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(3.14)       l'^*Gl'^u  and  L^*GS{T)y^  are  in  H^(E) 

2  9 
for  any  u  €  L   (E)   and         €  L  (Q). 

Next  we  show  that  L*FSi')Y^  and  L*FLu  are   in  H^(E)  for 

2  2 
y^  e  L  (Q)   and  u  €  L  (E).     In  fact,   for  y^  e  D(A),   we  have 

(3.15)  1^  FS(t)y^  =  FAS(t)y^  6  [0,T] ;  H^(Q)). 

By   (A. 5),   the  operator  FA  can  be  extended  to  be  bounded  on 

2  o 
L  (Q)    [1  p.   80].     Since  D(A)    is  dense  in  L   (Q),    the  formula 

2 

(3.15)   can  be  extended  to  all  of  L   (Q).     Consequently,  we 
deduce  that  v  - 

(3.16)  FS(.)yQ  e  H^(0)    for  y^  €  L^(Q). 

Let  f(s)   =  FS(s)y^  -  FS(T)y^  for  y^  €  L^(Q).  Then 

^1  2 
f  €  H   (0)    for  y^  €  L  (Q ) .     By  lemma  3.3   (v),  we  have 

(3.17)  L*f  e  hMe);    i.e.,   L*FS  (  •  )  yQ-L*FS  ( T )  y^  €  H^(E). 

On  the  other  hand,   by  theorem  1.2,   S* ( • -t ) FS ( T) y^  has  a 
well-defined  trace  on  r   in  H"^  (  [  t ,  T]  xr  )  .     In  fact,   by  lemma 
2.2,  we  obtain 

(3.18)  [D*(A*-kI)S*( .-t)FS(T)y^]  |^ 

=  a]^(  [S*(  .-t)FS(T)y^]  |^)~  which  is   in  H^([t,T]xr). 
From  (3.18),   for  0  <   t  <  T,  we  have 

(3.19)  (L*FS(T)y^) (t)    =   D*(A*-kI)    /J  S* ( s- t ) FS ( T ) y^ds 

rp 

=         A~   ( [S*(s-t)FS(T)y^] |^)"ds   is   in  H^(r). 
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Moreover , 

(L*FS(T)y^) (t)   =  D*  /J  S*(s-t) (A*-kI)FS(T)yQds 

=  D*  /^"^  S*(s) (A*-kI)FS(T)YQds 

is  absolutely  continuous  in  t.     Clearly,   the  derivative, 
-D*S* (T-t)A*FS(T)yQ  a.e.   t,   belongs  to  L^(Z).  Therefore, 
(3.19),  we  obtain  that  L*FS(T)y^  is  in  H-'-(E),  and 
consequently  from  (3.17)   that  L*FS{*)y^  is  in  H^(Z). 
Hence,   from  (3.14),  we  have    '      ' •  ' 

(3.20)  L*FS(.)y^  +  L^*GS(T)y^  6  H^(S)    for  y^  €  L^(Q) 

Next  we  will  show  that  L*FLu  €  ti^(Z)    for  u  €  L'^(Z). 
We  first  show  that  FLu  G       {Q)   for  u  €  L^(E).  Since 
FLU  e  C([0,T];   hJ^(Q)),   it   is  enough  to  show  that 
1^  FLu  e  L^(Q).     For  u   in  C^([0,T];  L^(r)), 

S( t-s)Du(s)ds  is  continuously  dif f erentiable  in  t 
[3,   p.   29].  Moreover,  ^        S ( t-s ) Du ( s ) ds  =  Du(t)  + 
A        S( t-s)Du( s)ds.     Since  F(A-kl)    is  bounded,  we  obtain 

(3.21)  ^  (FLu)(t)   =  F(A-kl)   ^         s(t-s)Du(s)ds  ^ 
=  F(A-kI)Du(t)   +  F(A-kl) (Lu) (t) 

+  F(A-kl)    k  /J  S(t-s)Du(s)ds   for  u  €  C^([0,T];  L^(r)). 

The  right  hand  side  of  (3.21)    is  in  L^(0)   for  u  6  L^(Z). 

Thus,   by  density  argument,  ~  (FLu)(t)   €  L^(0)  for 
u  e  L^iZ) . 
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That  is, 

(3.22)  Flu  e  H^{0)    for  u  €  L^(Z). 
Now  by  lemma  3.3   (v)   and   (3.22),  we  obtain 

(3.23)  L*(F.Lu-Fl'^u)   e  H^(Z)    for  u  €       ( Z  )  . 

By  the  similar  way  as  in  (3.17)  ~  (3.19),    it  is  easily 
verified  that   (3.24)   L*FLu  6  H^(S)    for  u  £  L^(S). 
Combining   (3.14)   with   (3.24),   we  obtain 

(3.25)        L*FL+l''^*Gl'^  maps  L^(S)    into  H-^(5:). 

Since  I  +  L*FL+l'^*Gl'^  has  an  inverse  on  L^(£),  from  (3.25)  we 
conclude  that  its  inverse  (  I  +  L*FL+l'^*GL^  )  maps  H^iZ)  into 
itself.     Thus,   by  (3.20),  we  obtain 

u°  =    (I  +  L*FL+l'^*GL^)~^  (L*FS(  .  )yQ  +  L^*GS(T)y^) 

is   in  H^(Z)    for  y^  €  L^(Q). 

This  completes  the  proof  of  the  theorem. 


CHAPTER  IV 
EVOLUTION  OPERATOR  AND  RICCATI  OPERATOR 

4.1     Evolution  Operator 

In  the  previous  chapter,   it  has  been  shown  that  the 

optimal  control  and  corresponding  trajectory  are  L^-functions 

2  2 

of  time  with  values   in  L  (r)   and  L  (Q)  respectively. 

However,   the  expression  for  the  optimal  control  in 
theorem  3.1   is  given  in  terras  of  the  initial  data.  This 
requires  knowledge  of  the  whole  process  starting  from  the 

....  u  .    .<  f  ... 

initial  time  t  =  0.  ' 

In  order  to  obtain  a  feedback  realization  of  the  optimal 
control  (pointwise  in  time),  we  need  a  representation  of  the 
optimal  control  as  an  operator  acting  on  the  trajectory  at 
each  time. 

In  order  to  accomplish  this,  we  introduce  an  evolution 
operator  $  which  describes  the  evolution  of  the  control  system. 
Let  s  be  an  arbitrary  time   in   [0,T).     We,  henceforth, 

regard  s  as  a  new  initial  time  for  our  control  problem,  with 

.    .  2 
new  initial  data  x  6  L  (Q). 

In  fact,   our  new  control  problem  over  the  time  interval 
[s,T]   with  initial  data  x  is  as  follows. 
(C.P)g     Minimize  the  functional  cost 
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2 

over  u  €  L  ([s,T]xr)   subject  to  the  mixed  problem 

(4.1)       a =  A( x,5 )y   in    [s,T] xQ 

My  =  u  in   [s  ,T]  xr 

y  (  s )  =  X  in  Q .  "     ,  ■  . 


Then,   by  the  same  argument  as  in  Chapter  III  for  (CP),  it 

is  clearly  seen  that  the  above  problem  (C.P)g  admits  a 

2 

unique  optimal  control   in  L  ([s,T]xr). 

Let  us  denote  by  u°(»,s;x)   and  y°(«,s;x)    the  optimal 
control  and  corresponding  trajectory  for  (C.P)g  respectively. 

In  order  to  obtain  explicit  representation  for 
u°(',s;x)   and  y°(»,s;x),  we  define  the  following  operators 
as   in  Chapter  II:  ,         •     ,  ■     . -.; 

(4.2)  (L  u)(t)   =   (A-kl)   jl  S( t-r)Du( r)dr  for  u  €  D(L  ), 

^  s  s 

(4.3)  (L*f)(t)   =     j    (L*f)(t),    s  <   t  <  T 

for  f  e  D( L* )  . 
0  ,0<t<s  ^ 

For  0  <   s  <   t  <  T,  we  define  the  operator  as 

L^u  =   (L^u)(t)    for  u  e  D(L^). 

^  s  s  I 

i 

Then,   by  the  same  argument  as   in  Chapter  II,   we  obtain  the 
following  lemma. 


50 


Lemma  4.1. 

(i)  The  operator         is  bounded   from  L'^([s,T]xr)  into 

C([s,T];    L'^(Q))   uniformly   in  s  €  [0,T]. 
2  t 

(ii)  For  each  u  €  L  (E),  L^u  is  continuous  in  both 
s  and  t,   0  <   s  <   t  <  T. 

(iii)  The  solution  yg  of   (4.1)    is  given  by 

y  (t)   =  S(t-s)x  -  L*^u  for  s  <   t  <  T. 
s  s 

Proof .     Part   (i)   and   (ii)    follow  directly  from  theorem 
2.2  and  corollary  2.2. 

For  part   (ii),    it   is,    in  view  of  theorem  2.2  (i), 

t  t  =;  2 

sufficient  to  show  that  L^u  =  L  u  -  S(t-s)   L  u  for  u  €  L  (S) 

2 

and  0  <   s  <   t  <  T.     In  fact,   for  u  €  L   (2),   0  <   s  <   t  <  T, 

lS  =   (A-kl)   jl  S(  t-r)Du(  r)dr 

=   (A-kl)         S( t-r)Du( r)dr  -   (A-kl)         S ( t-r ) Du ( r ) dr 
=  L^u  -   (A-kl)         S( t-r)Du(r)dr. 

We  next  show  that   (A-kl)         S ( t-r ) Du ( r ) dr  =  S(t-s)  L^u 
for  each  0  <   s  <   t  <  T.     For  x  in  D(A*),  we  have 

((A-kl)  S(t-r)Du(r)dr,x)^ 

=   (/^  S(s-r)Du(r)dr,   S*  ( t-s  )  (  A*-kI )  x ) 
=   {Jl  S(s-r)Du(r)dr,    ( A*-kI ) S * ( t-s ) x ) ^ 

=   (L^u,   S*(t-s)x)_   =   (S( t-s)L^u,x)^ 
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which  can  be  extended  to  all  of  L  (Q).     Thus,  we  obtain 
for  0  <   s  <  t, 

(A-kl)         S( t-r)Du( r)dr  =  S(t-s)L^u. 

This  completes  the  proof  of   the  lemma. 

Now  applying  the  procedures  of  Chapter  III  to  {C.P)g  , 
we  have 

(4.5)  u°(.,s;x)   =  L*Fy°(»,s;x)   +  L'^*Gy°(  T  ,  s  ;  x )  , 

s  s 

(4.6)  y°(.,s;x)   =  S(»-s)x  -  L^u° ( *  , s ; x) . 

We  obtain  the  following  lemma  similar  to  theorem  3.1. 

Lemma  4.2.     The  optimal  control  for  (C.P)^  is  given  by 
the  following  formula; 

u°(.,s;x)   =   (I+L*FLg+Lg*GLg)"^ [L*FS( .-s)x  +  Lg*GS ( T-s ) x] . 

From  (4.6),  we  are  in  a  position  to  define  an  evolution 
operator  as 

(4.7)  $(t,s)x  =  y°(t,s;x)  .        '  • 

for  0<s<t<T,xeL^(2), 

that   is,  $(t,s)x  =  S(t-s)x  -   ( LgU°( • , s ; x ) ) ( t ) . 

Then  $(t,s)    is  the  operator  which  describes  the  evolution  of 

the  optimal  trajectories. 

Some  properties  of  $  are  collected  below. 
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Lemma  4.3.     The  operator  <5  satisfies 

(i)  $(t,t)   =  I   identity  operator  for  0  <   t  <  T,  '" 

(ii)  •i(r,t)$(t,s)   =<I>(r,s)    for  0<s<t<r<T, 

(iii)  $(t,s)    is  strongly  continuous  in  t  from  L  (a)  .; 
into  C( [s,T] ;   L^(Q))    for  each  0  <   s  <  T, 

(iv)  ll'JCtjS)!!    „  <  c,   all  0  <   s  <   t  <  T 

L^(Q)-vL   (Q)  ■■.  ''  \ 

for  a  constant  c,  .■- ■ 

(v)  $(t,s)    is  strongly  continuous   in  s  from  L  (Q)  into 
C([0,t];   L^(Q))    for  each  0  <  t  <  T. 

Proof .     Part  (i)   is  obvious  from  (4.7).     Part  (ii)  can 

be  easily  obtained  from  the  uniqueness  of  optimal  control 

and  trajectory  for  a  given  data.     Part   (iii)    follows  from 

the  fact  that  the  semigroup  is  strongly  continuous  and  the 

operator  Lg  is  bounded  from  L^([s,T]xr)    into  C([s,T];  L'^(Q)). 

For  part   (iv),  we  first  show  that   |u  (.,s;x)L       '       <  c  Ixl 

'  I  [ s ,T] xr  I    I  a 

uniformly  in  s  for  all  x  €  L^(a).  :      ...  ..^ 

To  do  this,   we  use  lemma  4.2.     Since  L*FL  +L-'-*GL-'-   is  a 

s     s     s  s 

positive-definite  operator  on  L^([s,T]xr)   for  each  s,  we 
have  •  ?  ;• 

(4.8)  III  +  L*FL  +l'^*Gl'^II    ^  >  i 

ssss2  ^    J-  • 

L  (  [s,T]  xr)   -»>  itself 

Thus, 

(4.9)  II  (I  +  L*FLg+Lg*GLg)"^ll  <  1 

L  (  [s,T]  xr)  itself 

for  all  0<s<T[17].  • 
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Lemma  4.2  and   (4.9)  yield 

[s,T] xr 


(4.10)        |u°(. ,s;x) L  ^  < 


|L*FS(.-s)x  +  lT*GS(T-s)x|j^^^^^^ 

2 

for  X  e  L  (Q)   and  0  <   s  <  T. 

On  the  other  hand,    in  view  of  lemma  4.1   the  operator 
T 

norms  of  L*  and  L  *  are  uniformly  bounded  in  s;   i.e.,  for 

s  s 

0  <   s  <  T, 


(4.11)        1IL*II  <   ilL*ll    „  „ 

L^( [s,T]xQ)^L^( [s,T]xr)  L^(0)^L^(Z) 


L^(Q)^L^(  [s,T]xr)  L^(Q)^l'^(Z)  . 


By  (4.10)   and  (4.11),   together  with  the  boundedness  of 
FS('-s)   and  GS(T-s)   uniformly   in  s,   we  obtain 

(4.12)  'S'^M  [s,T]xr  ^  ^   I'^Iq  ^  ^  ^^^^^ 
and  0  <   s  <  T.     From  (4.2),   it  follows  that  for  0  <  s  <  T, 

(4.13)  IILgll  -  o  < 

L   (  [s,T]  xr)-^C(  [s,T]  ;   l'^(Q)  ) 

IlLII 

i  L   (S  )>C(  [0,T]  ;    l'^(Q)  )  . 

Thus,   from  (4.12)   and   (4.13),  we  have 

(4.14)  I  (L  u°(.  ,s;x)  )  (  t)  I      <   c  IxL 
for  0  <   s  <   t  <   T,   X  e  L^(Q)  , 
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where  c  is  a  constant  independent  on  s,   t  and  x.  Since 
|S(t-s)x|^  <   c   |x|^   for  0<s<t<T,   xeL(Q)   and  some 
constant  c,   from  (4.7)   and   (4.14)   we  obtain   (iv).     For  part 
(v),  we  recall  that  strong  continuity  (iii)    in  the  first 
variable  t  together  with  uniform  boundedness   (iv)  yields 
strong  continuity  in  the  second  variable  s  for  $(t,s)   as  in 
[1]  . 

This  completes  the  proof  of  the  lemma. 


4.2     Riccati  Operator  and  Feedback  Representation 

of  Optimal  Control 

In  this  section,  we  introduce  an  operator  (so  called  a 

2 

Riccati  operator)  on  L  (Q)  which  synthesize  the  optimal 
control   in  a  feedback  form.     Then  we  show  that   the  above 
operator  satisfies  a  Riccati  integral  equation. 

Let  us  define  an  operator  P(t)   on  L  (Q)   as  follows: 

(4.15)        P(t)x  =  /J  S*(s-t)F<5(s,t)xds  +  S*(T-t)G$(T,t)x 

for  xeD(P(t))    cL^(Q),   and  0  <   t  <  T. 

It  is  easily  seen,   from  lemma  4.2,   that  P(t)    is  a 

bounded  linear  operator  on  L  (Q)   for  each  t.     Moreover,  P(t) 

is  strongly  continuous  from  L^(Q)    into  C([0,T];   L^(Q))  with 

uniform  bound;    i.e.,   IIP(t)ll    ^  <   c  for  all  t  €  [0,T] 

L^(Q)-»-L^(Q) 

Now  we  are  in  a  position  to  prove  that  the  above 
operator  P(t)   satisfies  theorem  4.1. 


55 

Proof  o£  theorem  4.1   (i).     For  any  x  €  L'^(Q)  and 
0  <   s  <   t  <  T,   from  (4.5),   we  obtain 

(4.16)  u°(t,s;x)   =   (L*Fy°(« ,s;x) ) (t)   +   ( L^*Gy°( T , s ; x ) ) ( t ) 

=  D*(A*-kI)  S*(r-t)Fy°(r,3;x)dr 
+  D*(A*-kI)S*(T-t)Gy°(T,s;x) . 
Since  y°(r,s;x)   =  -I"  (  r ,  t  )$  ( t ,  s  )  x  ,    (4.16)  becomes 

(4.17)  u°(t,s;x)   =  D*(A*-kI)         S* ( r- t ) F$ ( r , t ) $ ( t , s ) xdr 

+  D*(A*-kI)S*(T-t)G$(T,t)$(t,s)x 
=  D*(A*-kI)P(t)$(t,s)x. 

Thus,   by  taking  t  =  s  and  x  =  y°(t)   =  $(t,0)y^,   we  obtain, 
from  (4.17),  that 

•      -  ■  ■ 

■  .    -  ■  .  i  .  .IS  ■ 

u°(t)    =  D*(A*-kI)P(t)y°(t)   =  BP(t)y°(t),  • 

where  B  =  D*(A*-kI)   is  given  in  lemma  2.2. 

This  completes  the  proof  of  theorem  4.1  (i). 

V'Je  first  prove  the  Riccati  integral  equation  in  theorem 

4.1   (ii)   under  an  additional  smoothness  assumptions  on  F  and 

G  in  this  section.     Then,   later  in  Chapter  V,  we  will  prove 

the  Riccati  equation  without  smoothness  assumptions  on  F  and 

G,   by  using  limit  process. 

Thus,    for  the  time  being,   let  us  assume 

2 

(A. 6)     F  and  G  map  L  (Q)   into  the  domain  of  A*. 
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Proof  of  theorem  4.1   (ii)   under  (A. 6).     We  recall,  from 
(4.7)   and  (4.17),   that  for  x  €  L^(Q)   and  U  <   s  <   t  <  T, 

(4.18)  <5(t,s)x  =  S(t-s)x  -   (  LgBP(  •  )$(  •  ,s)x)  (  t) 

=  S(t-s)x  -   (A-kl)         S( t-r)DBP(r)$( r,s)xdr. 

s 

2 

From  (4.15),  we  obtain  that  for  x,   y  €  L  (Q)   and  0  <   t  <  T, 

(4.19)  (P(t)x,y)^   =  (c5(s,t)x,FS(s-t)y)^ds 

+  (§(T,t)x,GS(T-t)y)^. 
Substituting  (4.18)   into  (4.19),   (4.19)  becomes 

(4.20)  (P(t)x,y)2  =  /J  ($(s,t)x,F-5(s,t)y)2ds 

+   ($(T,t)x,G$(T,t)y)^  ^   ^1  ^2 
where  ,  . 

1-1^  =         (F$(s,t)x,    (A-kl)   /J  S(s-r)DBP(r)$(r,t)ydr)^ds 

and  I2  =   (G$(T,t)x,    (A-kl)         S  ( T-r  )  DBP  (  r  )$  (  r ,  t )  ydr )  . 

Using  the  fact  that  (A*-kI)F  and  (A*-kI)G  are  bounded 

2 

on  L  (Q)   by   (A.6),   I-,^  and  1 2  become  . 

^1  "  -^t  ( (A*-kI)F$(s,t)x,  ll  S(s-r)DBP(r)$(r,t)ydr)^ds 

and  I2  =   ( {A*-kI)G$(T,t)x,         S ( T-r ) DBP ( r )$ ( r , t ) ydr ) ^ 
where  due  to  (A.6),  both  integrands  in        and  I2  are 
absolutely  integrable. 

Thus,  we  can  use  Fubini's  theorem  to  obtain         '■  [ 
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■    (4.21)         =  ( (A*-kI)F$(s,t)x,b(s-r)DBP(r)$(r,t)y)^drds 

=  /'J  (BS*(s-r)F$(s,t)x,BP(r)$(r,t)y)j,drds 

.     and  =  /J  ( (A*-kI)G$(T,t)x,S(T-r)DBP(r)$(r,t)y)2dr 

=  (BS*(T-r)G$(T,t)x,BP(r)$(r,t)y)j,dr. 
After  changing  the  order  of  integration  in  (4.21),   1-^  becomes 

■""1  "  ■'"t  ^1  (BS*(s-r)F$(s,t)x,BP(r)$(r,t)y)j,dsdr. 
Using  Fubini's  theorem  again,   we  obtain 

^1  "  -^t        •''r  S*(s-r)F$(s,r)$(r,t)xds,BP(r)$(r,t)y)j,dr. 
Thus,  we  have,   from  (4,15),  that 

(4.22)       ^1  +  ^2  "  •''t         ^1  S*(s-r)F$(s,r)$(r,t)xds 

+  BS*(T-r)G$(T,r)$(r,t)x,   BP ( r ) $ ( r , t ) y ) ^dr 
=  /J  (BP(r)$(r,t)x,BP(r)$(r,t)y)j,dr. 

From  (4.20)   and   (4.22),    it   is  readily  seen  that  P(t) 
satisfies  the  first  Riccati  integral  equation. 

This  completes  the  proof  of  theorem  4.1   (ii)  under  (A. 6). 

Notice  that  we  use  the  assumption   (A. 6)   only  to  justify 
change  of  the  order  of  integrations  in  l-j^  and 

Remark .  From  the  first  Riccati  integral  equation,  it 
follows  immediately  that  under  (A. 6),  our  Riccati  operator 
P(t)    is  self-adjoint  and  positive-definite. 


CHAPTER  V 
RICCATI  EQUATIONS 

5.1     Derivation  o£  Riccati  Equations 

The  Riccati   integral  equation  stated  in  theorem  4.1  is 
expressed  in  terms  of  the  feedback  evolution  operator  $.  " 
Another  representation  (so  called  the  second  Riccati 
integral  equation)  of  the  Riccati  operator  P(t)   in  terms  of 
only  the  original  semigroup  S(t)    is  available  under 
assumption  (A. 6).     Also,  a  Riccati  differential  equation  for 
P(t)    is  derived  under  assumption   (A. 6). 

The  purpose  of  this  section  is  the  proof  of  theorem 
5.1.     In  order  to  prove  the  theorem,  as  a  preliminary  step 
we  need  some  results  on  traces  of  the  Riccati  operator  P(t) 
as  well  as  on  regularity  of  the  evolution  operator  $. 

Remark.  The  presence  of  the  constant  k  in  the  operator 
B  does  not  affect  regularity  of  BP(t)  and  $ .  Thus,  from  now 
on,  without  loss  of  generality,   we  assume  k  =  0. 

Lemma  5.1.     Under  assumption   (A. 6), 

(i)  BP(.)x   is   in  C([0,T];   L^(r))    for  x  €  L^(Q), 

(ii)  BP(.)S(.)B*u   is  in  L^(Z)    for  u  6  L^(r). 

58 


59 


Proof  of  part  ( i ) .     From  the  assumption  (A. 6),  we 

2 

deduce  that  BS*(s-t)F  and  BS*(T-t)G  are  bounded  from  L  (Q) 
into  L^{r)  uniformly  in  s  and  t,  0  <  t  <  s  <  T.  Thus,  for 
X  €  L^(Q),   0  <   t  <  T,  we  have 

(5.1)    |BP(t)x|^  =    |/^  BS*(s-t)F$(s,t)xds+BS*(T-t)G$(T,t)x| 

<  /'^  |BS*{s-t)F<I.(s,t)x|^ds+|BS*(T-t)G$(T,t)x|^ 

<  /"J^  IIBS*(s-t)FII  ll'I>(s,t)  y  |x|^ds+IIBS*(T-t)Gll  ll$(T,t)  y  |x 

Combining   (5.1)  with  lemma  4.3   (iv),  we  obtain 

(5.2)  |BP(t)x|^  <   c   |x|^  uniformly  in  t  €  [0,T] 

for  x  e  L^(Q ) . 

The  continuity  of  BP(t)x  in  t  follows  from  strong  continuity 
of  S*  and  Thus,  part  (i)   is  proved. 

Proof  of  part  ( i  i )  .     We  recall,   from  (4.7)   and  lemma 
4.1,  that 

(5.3)  $(s,t)x  =  S(s-t)x-L^(I  +  LjFL^.  +  L^*GLj)"-^ 

[L*FS( •-t)x+L^*GS(T-t)x] 
for  xeL^(Q),0<t<s<T. 

We  next  show  that  b/'J  S*  (  s- t )  F$  (  s  ,  t )  S  ( t )  B*uds  and 

BS*(T-t)G$(T,t)S( t)B*u  are  well-defined   in  L^(E)  for 

2  0 
u  e  L  (D.     By  (5.3),  we  obtain  that  for  u  €  L  (r),   0  <   t  < 
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(5.4)       b/'J  S*(s-t)F$(s,t)S(  t)B*uds  =  b/^  S*  (  s- t )  FS  (  s  )  B*uds 

-   b/J  S*(s-t)FLf  (I  +  L*FL. +lI*Gl3)~''' [L*FS(  •  )B*u+L'^*GS(T)B*u]ds 
■'t  ttttt  t  t 

In  view  of   (A. 6)   and  boundedness  of  L*,   the  first  integral 

2 

on  the  right  hand  side  of   (5.4)    is   in  L  (Z).     From  (A. 6),  it 

is  easily  seen  that  FS(»)B*u  €  L^(0)   and  GS(T)B*u  6  L^(Q) 
2 

for  u  6  L  (r).     The  second  integral  of  the  right  hand  side 
of   (5.4)   becomes  '  ' 

D*/^  S*(  s-t)A*FL^(I  +  L*FLj.  +  L'J*GL^)~-^  [  L*FS  (  •  )  B* u+l'J*GS  (  T )  B*u]  ds 
2 

which  IS  in  L  (2)   since  the  integrand  is  uniformly  bounded 
in  t  e  [0,T].     Thus,   from  (5.4),  we  have 

b/'J  S*(s-t)F$(s,t)S(  t)B*uds  €  L^(E)    for  u  €  L^(r). 

Also,   from  (5.3),  we  have  that  for  u  €       ( r ) ,   0  <   t  <  T , 

(5.5)        BS*(T-t)G$(T,t)S( t)B*u  =   BS* ( T-t ) GS ( T ) B*u 

-  BS*(T-t)GL"J(I  +  L*FL,.+L^*GL^)"^  [LjFS(.  )B*u+L^*GS(T)B*u]  . 

By  the  same  argument  as  before,  we  obtain  that 

BS*(T-t)GS(T)B*u  €        ( S )    for  u  £  L^(r). 

The  second  part  of  the  right  hand  side  of  (5.5)   is  also  in 
2 

L  (E)   since   it  is  bounded  uniformly  in  t  €   [0,T].  Thus, 
from   (5.5),   BS* ( T-. ) G$ ( T , ♦ ) S ( • ) B*u  g  L^(E)    for  u  £  L^(r). 
Consequently,   we  conclude  that,   from  the  definition  of  P(t), 

BP(.)s(.)B*u  e  L^(z)  for  u  e  L^(r).  C 
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This  completes  the  proof  of  the  lemma. 
Vie  have  the  following  regularity  for  <5 . 

Lemma  5.2     Differentiability  of  $  in  [D(A*)]'. 
2 

For  X  e  L   (Q)   and  a.e.   t  6   [0,T],  we  have 

^<l>(s,t)x  =  -$(s,t)A[I-DBP(  t)  ]  X  where  the  equality  is 

understood   in  the  dual  space  [D(A*)]'. 

Proof .     We  first  show  that  for  x  €  L^(Q),  '. 
0  <   t  <   s  <  T, 

(5.6)       $(s,t)Ax  e   [D(A*)]'    uniformly  in  t  and  s.  > 

For  X  e  D(A),   0  <   t  <   s  <  T,   from  (4.7),  we  have 
(5.7)       $ ( s , t ) Ax  =  S ( s-t ) Ax  •  • 

-  L^(I  +  L*FLj.  +  L^*GL^)"^[L*FS(.-t)Ax  +  L^*GS  (  T- 1 )  Ax  ] 

=  AS(s-t)x-L^(I  +  L*FL^.+L^*GL^)"^  [  L  JFAS  (  • -t )  x+L^*GAS  (  T- 1 ) 

Since  FA  and  GA  can  be  extended  to  be  bounded  on  L^(Q)  in 

view  of   (A. 6),   L*FAS(--t)x+L^*GAS(T-t)x  £  L^([t,T]xr)  for 
2 

X  €  L  (Q).     Moreover,  we  have 

(5.8)        I (I+L*FL^+lJ*GL^)"^ [L*FAS( •-t)x 

+  Lj*GAS(T-t)x] I 

L^([t,T]xr) 

<   c   |x|     uniformly  in  t  €   [U,T]    for  x  €  L^(Q) 
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On  the  other  hand,   since         is  bounded  from 
L^([t,T]xr)    into  C([t,T];   L'^(Q)),   we  deduce  from  (5.8)  that 

(5.9)  the  second  terra  in  the  right  hand  side  of  (5.7) 

2 

is   in  C([t,T];   L   (Q))   as  a  function  of  s, 
uniformly  in  t  6  [0,T]. 
It  is  clear  that 

(5.10)  S(»-t)Ax  e  C([t,T];   L^(Q))    for  x  £  D(A),  _ 
S(»-t)Ax  e  C([t,T];    [D(A*)]')    for  x  6  L^(Q), 


uniformly  in  t  €  [0,T]. 

Hence,   from  (5.7),    (5.9)   and   (5.10),  we  conclude  that 

(5.11)  $(»,t)Ax  €  C([t,T];   L^(Q))    for  x  €  D(A), 
$(»,t)Ax  e  C([t,T];    [D(A*)]')    for  x  €  L^(Q) 

uniformly  in  t  €   [0,T].  .  . 

/ 

From  (5.11)   and  lemma  5.1   (i),  we  arrive  at 

(5.12)  $(• ,t)ADBP(t)x  e  C( [t,T] ;  [D(A*)]') 

2 

for  x  6  L  (Q),   uniformly  in  t  €  [0,T]. 

with  the  aid  of  (5.12),  we  now  justify  the  following 
calculations   in  the  space   [D(A*)]'   as  in[l]. 
For  X  e  L^(Q)   and  a.e.    t  €    [0,T],    0  <   r  <   s  <  T, 

(5.13)  0  =  ^  $(s,r)x  =  1^  [$(s,t)$(t,r)x] 
=  $(s,t)]$(t,r)x  +  $(s,t)  l^$(t,r)x. 
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Since  ^(t,r)   is  the  feedback  evolution  operator, 

^$(t,r)x  =  A'I>(t,r)x-ADBP(t)$(t,r)x       a.e.   t  €  [r,T]. 

Thus,    in  view  of   (5.11)   and   (5.12),    (5.13)  becomes 

(5.14)  0  =  ^(s,t)]$(t,r)x 

+  $(s,t) [A$(t,r)x-ADBP(t)$(t,r)x]    in  [D(A*)] 
Hence,  we  have 

(5.15)  $(s,t)]$(t,r)x  =  -<$(s,t)A[I-DBP(t)  ]$(t,r 

for  X  e  L^(Q)   and  a.e.   t  €   [0,T] ,   r  <   t  <  s,  where 
by  (5.12),   the  right  hand  side  is  well  defined  in  [D(A*)] 
Thus,   the  lemma  follows  from  (5.15)   after  taking  r  = 

Lemma  5.3.     For  x  €  L^(Q)   and  0  <   t  <  s  <  T, 

$(s,t)x  =  S(s-t)x  -  $(s,r)ADBP(r)S(r-t)xdr. 
Proof .     From  lemma  5.2,  we  have 

(5.16)  $(s,t) ] S( t-r)x  e   [D(A*)]'    for  x  €  L^(a) 
and  a.e.   t  6   [0,T],   0<r<t<s<T.     That   is,  <' 

(5.17)  -$(s,t)AtI-DBP(t)]S(t-r)x  €  [D(A*)]' 
2 

for  X  e  L  (Q)   and  r  <   t  <   s  uniformly  in  t  €  [0,T]. 

Since   (5.17)    is  uniformly  bounded   in  t  6   [0,T]  and 
[D(A*)]'    is  complete,   we  obtain  that 
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(5.18)  $(3,t)A[I-DBP(t)]S(t-r)xdt   is  well  defined 

in   [D(A*)]'    for  x  €  L^(Q)   and  0  <   r  <         <   t2  <   s  <  T . 
2 

That  IS,   for  x  6  L  (Q)   and  0  <   r  <   t ^   <         <   s  <  T , 

U  ^(s,t)]S(t-r)xdt  €  [D(A*)]'. 

On  the  other  hand,  after  integrating  by  parts,  we  obtain 
t  , 

(5.19)  /^.^  $(s,t)]S(t-r)xdt 

t  —  t  ^  t  ^ 

=   [$(s,t)S(t-r)]     /-  -  j       $(s,t)AS(t-r)xdt  ■ 

^1        ^1  :. 

where  integrations  are  performed  in  the  space  [D(A*)]'. 
From  (5.18),    (5.19)  becomes 

(5.20)  $(S,t^)S(tj^-r)x  =  $(  s,t2)S(  t2-r)x 

-  /  $(s,t)ADBP(t)S(t-r)xdt 
" 

2  ' 
for  x  e  L  (Q)   and  0  <   r  <   t^  <   t2  <   s  <  T. 

2 

By  taking  t^  =  r  and  t2  =  s,     we  obtain,   for  x  €  L  (Q), 

(5.21)  $(s,r)x  =  S(s-r)x  -        $ ( s , t ) ADBP ( t ) S ( t-r ) xdt . 

2 

Since  '5(s,r)x  and  S(s-r)x  are  in  L  (Q),   the  second  term  in 

2 

the  right  hand  side  of   (5.21)  must  be  in  L  (Q). 
This  completes  the  proof  of  the  lemma. 
Now  we  are   in  a  position  to  prove  theorem  5.1. 
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Proof  of  theorem  5.1(i).     For  x,   y  €  L^(Q)  and 
0  <   t  <  T,   from  (4.15)   we  have 

(5.22)  (P(t)x,Y)^   =  /J  ($(s,t)x,FS(s-t)y)^ds 

+  ($(T,t)x,GS(T-t)y)^. 
By  lemma  5.3,    (5.22)  becomes 

(5.23)  (P(t)x,y)^  =  (S(s-t)x,FS(s-t)y)ds 

+   (S(T-t)x,GS(T-t)y )^  "^1~^2 

where         =         (/^  $ ( s , r ) ADBP ( r ) S ( r- t ) xdr , FS ( s-t ) y ) ^ds 

and  I2  =   (/^  $(T,r)ADBP(r)S(r-t)xdr,GS(T-t)y)^. 

In  order  to  change  the  order  of   integrations   in  1-^  and 

2 

1-2,   we  need  to  show  that  for  x,   y  6  L  (Q), 

(5.24)  ($(s,r)ADBP(r)S(r-t)x,FS(s-t)y)^drds 
and  ,  ■ 

(5.25)  j1  ($(T,r)ADBP(r)S(r-t)x,GS(T-t)y)-dr 

are  Lebesque   integrals.     From  (5.12),    for  x  €  L^(Q)  we 
obtain 

(5.26)  $(s,r)ADBP(r)S(r-t)x  €  [D(A*)]' 
uniformly  bounded  in  r,   s  and  t. 

Thus,   from  (5.12)   and   (A. 6),    it  follows  that 

(5.27)  ($(s,r)ADBP(r)S(r-t)x,FS(s-t)y)^ 
and  ($(T,r)ADBP(r)S(r-t)x,GS(T-t)y) 
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are  well-defined  and  uniformly  bounded  in  r,  s  and  t. 
This   implies  that   (5.24)   and   (5.25)   are  Lebesque  integrals 
Now  we  are  in  a  position  to  apply  Fubini's  theorem   [6]  to 
and  Then  we  have 

...    (5.28)         =  /"J  /^($(s,r)ADBP(r)S(r-t)x,FS(s-t)y)^drds 
=  /J  /^($(s,r)ADBP(r)S(r-t)x,FS(s-r)S(r-t)y)^dsdr 
=  !l  j1  (S*(s-r)F$(s,r)ADBP(r)S(r-t)x,S(r-t)y)^dsdr, 


(5.29)         =  ('5(T,r)ADBP(r)S(r-t)x,GS(T-t)y)^^dr 


=  /J  (S*(T-r)G$(T,r)ADBP(r)S(r-t)x,S(r-t)y)^^dr. 


Similarly,   from  (5.12)   and   (A. 6),   it  is  easily  seen 

that 

/j,  S*  (  s-r )  F$  (  s  ,  r )  ADBP(  r  )S  (  r-t )  xds 
is  uniformly  bounded  in  r.     '  ' 
Hence,  by  Fubini's  theorem,   from  (5.28)  we  obtain 

(5.30)  =         (/^  S*(s-r)F$(s,r)ADBP(r)S(r-t)xds, 

S( r-t)y )Qdr. 

Combining  (5.29)  with  (5.30),  together  with  the  remark  in 
section  4.2,  we  have  that 


(5.31)  +       "  -^t  (P(j^)ADBP(r)S(r-t)x,S(r-t)y)^dr 

=  /"^   (BP(r)S(r-t)x,BP(r)S(r-t)y)  dr 
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which  is  well  defined  in  view  of  lemma  5.1   (i).     Thus,  the 
first  part  of  theorem  5.1   follows  from  (5.23)   and  (5.31). 

Proof  of  theorem  5.1    (ii).     The  Riccati  differential 
equation  is  obtained  by  differentiating  the  Riccati  integral 
equation  with  respect  to  t  for  x,   y  in  D(A).     In  fact,  for 
X,  y  in  D(A)   and  a.e.   t  in   [0,T],   by  differentiating  the 
second  Riccati   integral  equation  with  repspect  to  t,  we 
obtain 

(5.32)  ^  (P(t)x,y)^  =  -(x,Fy)2-/J  ( S ( s-t ) Ax , FS ( s- t ) y ) ^ds 

rp 

-  (S(s-t)x,FS(s-t)Ay)2ds+(BP(t)x,BP(t)y)^ 
+  (BP(s)S(s-t)Ax,BP{s)S(s-t)y)j,ds 

rp 

+  /  (BP(s)S(s-t)x,BP(s)S(s-t)Ay)„ds 

^  r 

-  (S(T-t)Ax,GS(T-t)y)^  -   ( S ( T- t ) x , FS ( T- t ) Ay ) ^ . 

All  terms  appearing  in  (5.32)   are  well  defined  in  L^-sense, 
in  view  of  lemma  5.1  (i). 

After  applying  the  second  Riccati  integral  equation  to  the 
inner  products   (P(t)Ax,y)^     and  (P(t)x,Ay)^  in  (5.32),   it  is 
easily  seen  that  for  x,  y  €  D(A)   and  a.e.   t  6  [0,T], 

^  (P(t)x,y)^   =  _(x,Fy)^  -   (P(t)Ax,y)^  -  (P(t)x,Ay)^ 
+   (BP(t)x,BP(t)y)^  with  P(T)   =  G. 

This  completes  the  proof  of  theorem  5.1. 
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5.2     Uniqueness  of  Riccati  Operator 

In  this  section,  we  give  the  proof  of  theorem  5.2. 
Let    V   be  the  class  of  self-adjoint,  positive-definite  and 
bounded  linear  operators  P(t),  0  <  t  <  T  on  L^(Q)   having  the 
properties  that 

(5.33)       BP(»)    is  bounded  from  L   (Q)    into  L° { [0 ,T] ;   L  (V) 

and  (5.34)       BP(.)S(»)B*  is  bounded  from  L^(r)    into  L^(Z). 

VJe  show  that  the  Riccati  equations  in  theorem  5.1  admit 
the  unique  solution  in    V.     Since  the  second  Riccati 
integral  equation  can  be  derived  from  the  Riccati  differ- 
ential equation  in  a  usual  way,   the  uniqueness  for  the 
Riccati  differential  equation  follows  from  the  uniqueness 
for  the  Riccati  integral  equation. 

Thus,   it  is  enough  to  prove  the  uniqueness  for  the  second 
Riccati  integral  equation  in  the  class    P    (theorem  5.2). 

Proof  of  theorem  5.2.     Since  our  Riccati  operator  P(t) 
defined  in  (4.15)   is  a  solution  to  the  Riccati  equations  in 
theorem  5.1  and  it  satisfies  (5.33)   and   (5.34),  the 
existence  of  a  solution  in    P    is  guaranteed. 

Suppose  that  an  operator  P(t)  e    V    is  also  a  solution 
to  the  second  Riccati  integral  equation. 
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In  fact,   for  x,   y  e  L  (Q)   and  0  <   t  <  T, 

(5.35)  (P(t)x,y)^   =  (S(s-t)x,FS(s-t)y)^ds 

T  __   

-  (BP(s)S(s-t)x,BP(s)S(s-t)y)j,ds+(S(T-t)x,GS(T-t)y) 

We  will  show  that  P(t)   =  P(t)   for  0  <  t  <  T. 
Let  us  consider  the  following  mixed  problem:  for 
y^  €  L^(Q) , 

(5.36)  d  ^y  =  A( x,5  )y       in  Q 

My  =  BP'(  •  )y  ( •  )  in  Z 
y(0)     =  y^  in  Q 

We  first  show  that  (5.36)   has  a  unique  solution  in  L^(0). 
This  can  be  done  by  Rauch's  estimate  (1.9)   and  a  fixed  point 
theorem.     VJe  rewrite   (5.36),  using  semigroup  representation, 

(5.37)  y  =  S(.)yQ  -  LBP(.)y(0. 

Let  us  define  an  operator  E  as  follows: 

(5.38)  (Ey)(t)   =   (LBP( . )y( • ) ) ( t) ,   0   <   t  <  T,  ' 

for  y  e  D(E)    c  l^(0) . 

Then  it  is  easily  seen,   from  theorem  2.2   (i)   and  (5.33), 
that  E  is  a  bounded  linear  operator  on  L^(Q). 
By  (5.38),    (5.37)  becomes 

(5.39)  y  =  s( .  )yQ  -  Ey. 
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In  order  to  show  that   (5.39)   has  a  unique  solution  in 
L°°([0,T];   L^(Q)),   we  show  that  E   is  a  contraction  map  on 
l"([0,T];   L^(Q))    for  small  T  >   0.     For  y  €  L°°{[0,T];  L^(Q)), 

—  00  2 

we  have,   from  (5.33),   that  BP(')y(')   €  L  ([0,T];   L  (Q)). 
Thus,   by  theorem  2.2   (i),    it  follows  that 

r-:       (5.40)        lEyl  „  <   c   I  BP(  •  )y(  •  )  I  ^ 

l'( [0,T] ;   L^(Q) )  L^(Z) 

^  <   c '   /T   I BP( • )y ( • ) I  ^ 

l"'(  [0,T]  ;    L^(r)  ) 

for  y  e  l"([0,T];  L  (Q))  where  c  and  c'   are  constants 
independent  on  sufficiently  small  time  T. 
By  (5.33),    (5.40)  becomes 

lEyl  <  c  /f  lyl  „ 

l"([0,T];    L^(Q))  l"([0,T];  L'^(Q) 

for  y  e  L°°([0,T];   L  (Q))  where  c  is  a  constant  independent 

on  small  T.     Thus,   E  can  be  made  a  contraction  map  by  taking 

Tq  >  0  sufficiently  small,   i.e.,  c  /T^  <  1. 

By  a  fixed  point  theorem,    (5.37)   has  a  unique  soultion 

y  in  l"([0,Tq];   L^(Q))    for  small  T^  >  0. 

Since  our  system  (5.36)    is  time-invariant  and  C-*- ,  it 
can  be  shown  that  for  any  given  T  >  0,    (5.36)   has  a  unique 
solution  by  repeating  the  above  argument  n  times  where  n  is 
an  integer  such  that  n  T^  >  T.     Let  u  =  BP(»)y(*). 

Our  next  step  is  to  show  that  "{J  = 

To  this  end,  we  introduce  the  evolution  operator  $ 
corresponding  to  (5.36)   as  in  Chapter  IV. 
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Then,  as  in  Chapter  IV,   from  (5.37),  we  have 

(5.41)  ¥(s,t)x  =  S(s-t)x-(L^BP( . )$( • ,t)x) (s) 

2 

for  X  €  L  (Q)   and  0  <   t  <   s  <  T. 

2 

That  IS,   for  x  €  L  (Q)   and  0  <   t  <   s  <  T, 

i'(s,t)x  =  S(s-t)x-A/^  S(s-r)DBP(r)$(r,t)xdr. 
Substituting   (5.41)    into  (5.35),  we  obtain  for  x,  y  e  L^(Q), 

(5.42)  (P(t)x,y)^   =  (¥(s,t)x,FS(s-t)y)Qds 

T  s 

+  /j.  (A/^S(s-r)DBP(r)$(r,t)xdr,FS(s-t)y)^ds 

-  /J  (BP(s)$(s,t)x,BP(s)S(s-t)y)j,ds 
•T   ,-,    .    .8  ^ 


-         (BP(s)A/^  S(s-r)DBP(r)$(r,t)xdr, 


BP(s)S(s-t)y)   ds  +  ($(T,t)x,GS(T-t)y) 


Q 


Q 


+  (a/^  S(T-r)DBP(r)¥(r,t)xdr,GS(T-t)y) 
=  !^  (S*(s-t)F$(s,t)x,y)^ds 

+   (S*(T-t)G¥(T,t)x,y)^  + 

T  s 

where         =         (a/^  S ( s-r ) DBP ( r )¥(  r , t ) xdr , FS ( s- t ) y ) ^ds , 

T       

''^2  "  -^t  (BP(s)$(s,t)x,BP(s)S(s-t)y)^ds, 

T      _  s  _  _ 

I3  =         (BP(s)a/j.  S(s-r)DBP(r)$(r,t)xdr,BP(s)S(s-t)y)j,ds 

and        =   (a/^  S(T-r)DBP(r)$(r,t)xdr,GS(T-t)y)Q. 
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Changing  the  order  of   integrations   (this  is  justified  in 
view  of  (A.6)),  we  obtain 

(5.43)  =  (S(s-r)B*BP(r)$(r,t)x,FS(s-t)y)2drds 

=  /t  /r  (S(s-r)B*BP(r)¥(r,t)x,FS(s-t)y)Qdsdr 

and   (5.44)  =         ( S (T-r ) B*BP ( r )¥( r-t ) x ,GS ( T-t ) y ) ^dr . 

Since  (BP(s)$(s,t)x,BP(s)S(s-t)y)p  =   ( P ( s ) B*BP ( s )¥ ( s , t ) x , 
S(s-t)y)2   ,  we  obtain 

(5.45)       I2  =  (P(s)B*BP(s)$(s,t)x,S(s-t)y)Qds. 

On  the  other  hand,   in  order  to  use  Fubini's  theorem  in  I3, 
we  need  to  show  that 

"^t  •''r  (BP(s)S(s-r)B*BP(r)$(r,t)x,BP(s)S(s-t)y)j,dsdr 

is  Lebesque  integral.     This  can  be  done  as  follows. 
In  fact,   by   (5.33)   and   (5.34)    together  with  Schwarz ' s 
inequality,  we  obtain  the  following  inequalities; 

/J!   I  (BP(s)S(s-r)B*BP(r)¥(r,t)x,BP(s)S(s-t)y)j,|dsdr 

1 

<  /J   [/^  |BP(s)S(s-r)B*BP(r)$(r,t)x|^ds] 2 

1^ 

t/^   |BP(s)S(s-t)y l^ds] ^dr 
T        I    — 

<  c   |BP(r)$(r,t)x|j,  |y  l^dr  for  some  constant  c  >  0 

<  c'    Ixl    lyl     for  some  constant  c'   >  0. 
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Thus,   by  Fubini's  theorem,  we  have 

(5.46)  I3  =  !^  il  (BP{s)S(s-r)B*BP(r)¥(r,t)x, 
BP( s)S( s-t )y ) pdrds 

=  /J  (BP(s)S(s-r)B*BP(r)¥(r,t)x,BP(s)S(s-t)y)j,dsdr. 

Let  u(s)   =  B*BP(s)$(s,t)x  and  v(s)   =  S(s-t)y. 
Then  u(.)   6  L"([t,T];    [D(A*)]')   and  v(.)   €  C([t,T];  L^(Q)). 
Vvfe  note  that  all  integrals  in  (5.43)  ~  (5.46)   are  Lebesque 
integrals.     That  is, 

^1  "  -^t  ^r  (S(s-r)u(r)  ,FS(s-r)v(r)  )gdsdr, 

12  =  !l  (P(r)u(r) ,v(r) )^dr, 

13  =        il  (BP(s)S(s-r)u(r)  ,BP(s)S(s-r)v(r)  )j,dsdr 

and        =  /J  (S(T-r)u(r) ,GS(T-r)v(r) )Qdr 

are  all  Lebesque  integrals. 

Combining  I^,   1^,         and         ,  we  obtain 

(5.47)  1^-12-13+1^=         [/^(S(s-r)u(r) ,FS(s-r)v(r) )Qds 

-  (P(r)u(r) ,v(r) )2-/ J( BP ( s ) S ( s-r ) u ( r ) , BP( s ) S ( s-r ) v ( r ) )^ds 
+(S(T-r)u(r) ,GS(T-r)v(r) )g]dr. 

On  the  other  hand,   in  view  of  (5.33)   and  (5.34),  the 

second  Riccati   integral  equation  holds  for  x  e   [D(A*)]'  and 
2 

y  e  L  (Q),   thus,   in  particular  for  u(r)   and  v(r),  a.e. 
re   [t,T] . 
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Applying  the  second  Riccati  integral  equation  to  (5.47)  with 
X  =  u(r)   and  y  =  v(r)   a.e.   r  €   [t,T],  we  obtain 


(5.48)  ^1  "  ^2  ~  -"-S  ^  ■'■4  "  °* 

From  (5.42)   and   (5.48),    it   is  readily  seen  that 

(5.49)  (P(t)x,y)^  =  (S*(s-t)F$(s,t)x,y)^ds 

+   (S*(T-t)Gi"(T,t)x,y)^   for  all  x,  y  in  L^(Q) 
and  t  in   [0,T].     Thus,   from  (5.49),  we  obtain 


(5.50)        P(t)x  =  /"J  S*(s-t)F$(s,t)xds  +  S*  (  T- 1 )  G$  ( T  ,  t ) 

2 

for  all  X  in  L   (Q)   and  t  in  [0,T]. 


Since  $  is  the  evolution  operator  for   (5.36),   for  t  <  s, 
y(s)   =$(s,t)y(t).     Therefore,    (5.50)  becomes 

(5.51)  P(t)7(t)   =  /"J  S*(s-t)F$(s,t)7(t)ds 

.  .  :        +  S*(T-t)G$(T,t)y(  t) 

=         S*(s-t)Fy(s)ds  +  S*(T-t)Gy(T) . 

From  (5.51),  we  obtain 

(5.52)  u(t)   =  BP(t)y(t)   =  b/^  S* ( s-t ) F$ ( s , t ) y ( t ) ds 
+  BS*(T-t)G$(T,t)7(t) 

=   (L*F7)(t)   +   (L^*Gy(T) ) (t)   for  a.e.   t  €  [0,T] 

From  (3.7),  (5.52)  implies  that  u  satisfies  the  optimality 
condition  for  our  control  problem. 
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By  the  uniqueness  of  the  optimal  control  to  our  control 
problem,  we  conclude  that  u  =  u*^. 

Thus,  y  =  y°  and  consequently,  $  =  $.     By   (4.15)   and  (5.50), 
we  conclude  that  P(t)   =  P(t)   for  all  t  €  [0,T]. 
This  completes  the  proof  of  theorem  5.2. 

Remark .     VJe  use   (5.33)   and   (5.34)   to  justify  change  of 
the  order  of  integrations  in  I^,   l^r   I3  and  I^,  and  also  to 
verify  the  second  Riccati   integral  equations  for  u(r)  and 
v( r)   a.e.  r. 

5.3  Regularization 

In  this  section,  we  assume  that  F  and  G  are  arbitrary 

self-adjoint,  positive-definite  bounded  operators  on 
2 

L  (Q ) ;   i.e.,  we  do  not  assume  any  smoothness  on  F  and  G. 

Vve  will  show  that  the  corresponding  Riccati  operator 

defined  in  (4.15)   can  be  determined  as  the  strong  limit  on 
2 

L  (Q)  of  a  sequence  of  solutions  to  regularized  Riccati 

equations  in  theorem  5.1.     In  fact,   the  purpose  of  this 

section  is  to  give  the  proof  of  theorem  5.3. 

Suppose  that  { F^}  ^  and  {g^}^  are  sequences  of  self- 

2 

adjoint,  positive-definite  operators  on  L  (Q)  which  satisfy 
(5.53)       F^  and  G^  map  L^(Q)    into  D(A*)    for  all  n, 

2 

F     ->•  F  and  G     -»■  G  strongly  in  L  (Q)   as  n  ». 
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Such  sequences  can  be  easily  constructed  as  we  will  see 
later.     For  simplicity  of  notation,   let  us  denote  by  u°, 
y^,  P^,  etc.,   the  optimal  control,   trajectory,  evolution 

operator,  Riccati  operator,  etc.,  corresponding  to  the 
smoothing  operators         and  G^,  in  the  functional  cost  (1.12). 

In  order  to  prove  theorem  5.3,  we  need  the  following 
lemma . 


2 

Lemma  5.4.     For  any  given  x  €  L  (Q), 

(i)  u°(»,t;x)  >  u°(»,t;x)    in  L^([t,T]xr)   as  n 

uniformly  in  t,   0  <   t  <  T, 

2 

(ii)  $^(s,t)x  +  $(s,t)x  in  L  (Q)   as  n  -»•  « 
uniformly  in  both  s  and  t,   0  <  t  <  s  <  T. 


Proof  of  part  ( i ) .     From  lemma  4.2,  we  recall  that  for 
each  n  and  0  <   t  <  T, 

(5.54)        u°(.,t;x)    =    [H-LJF^L^+lJ*G^L^]~^ [LjF^S( .-t)x 


+  l1*G  S(T-t)x] . 
t  n 


Let         =  I+LfF^L^+Lj*G^Lj  and  R^  =   I+L*FL^+L^*GL^ . 

-1-1  2 
We  will  show  that  R^     >  R^     strongly  in  L  ([t,T]xr)   as  n  >  » 

uniformly  in  t  €  [0,T]. 

2 

For  a  given  u  €  L  ([t,T]xr),     we  have 
(5.55)       IV-Vl[t,T)xr  <  l'-f<fn^t"-''V>l[t,Tlxr 

*  i^I*<=n^X"'l[t,Tlxr. 
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In  view  of  lemma  4.1   (i),   (ii),   it  is  easily  seen  that 
T 

L*  and  L^*  are  uniformly  bounded  in  t.     Thus,    (5.55)  becomes 

(5.56)       |r  u-R  u|  <  c   If  L^u-FL^uL^  ^1  ^ 

In       o   l[t,T]xr  I    n  t         t  i[t,T]xQ 

I        T         T  I 
+  c    G  L^u-GL^uK    for  some  constant  c  >  0. 
I    n  t         t   I Q 


We  recall  that  strong  convergences  of  bounded  linear 

operators  on  a  Hilbert  space  turn  to  be  uniform  convergences 

on  compact  subsets  of  the  space.     Thus,   from  (5.53),   F^  >  F 

2 

and  G|^  ->■  G  uniformly  on  compact  subsets  of  L  (Q)   as  n  =>. 

On  the  other  hand,   from  lemma  4.1   (ii),  we  deduce  that 

{l^u   I   0  <   t  <t}   and  (l^u   |  0  <  t  <  s  <  t}   are  compact 
2 

subsets  of  L  (Q),   being  the  images  of  compact  sets  under 

continuous  maps.     Hence,  we  obtain 

(5.57)    I^j^l'^u-GL^uI^  >  0  as  n  ->  00  uniformly  in  t  €  [0,T] 


and 


I  Fj^L^u-FL^u  1 2  ->•  0  as  n  ■>  »  uniformly  in 
both  s  and  t,0<t<s<T. 


Combining  (5.56)  with  (5.57),  we  conclude 


(5.58)        i^n'^~^o"l  [t,T]  xr       0         n      «  uniformly  in  t; 

2 

I.e.,  ^  strongly  in  L  ([t,T]xr)  as  n  ->  o  uniformly  in 
For  u  €  L^( [t,T] xT ) , 

lC-^;'-l[t,t]xr  =  lC(V^o)C-l[t,T]xr 

I  ^^n'^o^^o"*^^!  [t,T]  xr       0         n  +  »  uniformly  in  t, 
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where  in  the  last  inequality,  we  use  the  fact  that  H  Rj^"*" "  <  1 
for  all  n,  and  (5.58).     Thus,  we  obtain 

(5.59)       r""""  >  r""^  strongly  in  L^([t,T]xr)   as  n  •>  « 

uniformly  in  t  €  [0,T]. 


Let  V     =  L*F  S(»-t)x  +  lI*G  S(T-t)x, 

n        t  n  t  n 


V  =  L*FS(»-t)x  +  L'J^*GS(T-t )  X. 


Then  in  a  similar  manner,   it  is  easily  seen  that 

2 

(5.60)  >  V  in  L  ([t,T]xr)  as  n  -»■  »  uniformly  in  t. 
From  (5.54),  we  obtain 

(5.61)  |u^(.,t;x)-u°(.,t;x)|^^^^j^^  =    I  C^n'^^'H  [  t  ,T]  xF 

'   l^'vCH[t,T]xr  ^  l^n'^-^o'H[t,T]xr 


<    IIR~"'"ll|v  -vL^         „  +    I  R~''-v-R~''-v         ^,  „. 
n     In     I [t,T] xr       In         o     I [t,T] xr 


Since  IIR~^II    <   1   for  all  n,   by  (5.59)   and   (5.60),   the  right 

hand  side  of  (5.61)  converges  to  0  as  n  •>  ®  uniformly  in 
t  €   [0,T] . 

This  completes  the  proof  of  part  (i). 


Proof  of  Part  ( ii ) .     For  t  <   s  and  x  €  L   (Q),  from 
(4.7)   we  obtain 
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(5.62)  $^(s,t)x-$(s,t)x  =   (L^(u°( . ,t;x)-u°( . ,t;x) ) ) (s) 

2  9 
Since  L,.   is  bounded  from  L   ([t,T]xr)    into  C([t,T];   L  (Q)) 

uniformly  in  t,    (5.62)  yields 

(5.63)  |*^(s,t)x-$(s,t)x|^ 

<  |l^(u°(. ,t;x)-u°(. ,t;x) ) I 

C( [t,T] ;   L^(2) ) 

<  c    |u°( . ,t;x)-u°( . ,t;x) | 

L   (  [t,T]  xD 

for  some  constant  c  >  0  where  the  last  term  converges  to  0 
as  n  ^  o»  uniformly  in  t  by  part  (i). 

This  completes  the  proof  of  part  (ii). 

Remark .     From  lemma  5.4,  we  infer  that 

(i)  u°  *  u°  in  L^(r)  as  n  >  », 

(ii)  y^^t)  -»■  y  (t)  uniformly  in  t  as  n  ->■  00. 
Now  we  are  in  a  position  to  prove  theorem  5.3. 

Proof  of  theorem  5.3.     For  any  t  €   [0,T]  and 
X  e  L^(Q ) ,     we  have 

(5.64)        |P^(t)x-P(t)x|^  I 

<   Jl   |s*(s-t) [F^$^(s,t)x-F$(s,t)x] l^ds 
+   |s*(T-t)  [G^$^(T,t)x-G<5(T,t)x]  1^. 
Since  S*(t)   is  a  strongly  continuous  semigroup,   its  norm  is 
bounded  on  a  finite  time   interval.     Thus,    (5.64)  becomes 


80 


(5.65)        |p  (t)x-P(t)x|     <   c/I"   |f$   (s,t)x-F$(s,t)x|  ds 

+  c    |G^$^(T,t)x-G5  (T,  t)x|^   ^       i't  l^n^n^^'^^^'^n^^^'^^^lQ^^ 
+  c         I F^$(s,t)x-F$(s,t)x|^ds  +  c  |G^$^(T,t)x-G^$(T,t)x|Q 
+  c    |g  $(T,t)x-G$(T,t)xl 

for  some  constant  c  >  0.     From  lemma  4.3   (iii)   and   (v),  we 
deduce  that 

(5.66)  {$(s,t)x   I    0  <   t  <   s  <  t}  and 

{$(T,t)x   I   0  <  t  <  t}   are  compact  subsets  of  L^(Q) 

Thus,   by  the  similar  argument  to  the  proof  of   (i)    in  lemma 
5.4,  we  obtain  that  by  (5.66), 

(5.67)  |F^$(s,t)x-F$(s,t)x|^  0 
uniformly  in  both  s  and  t  as  n  ^  », 

and  |g  $(T,t)x-G$(T,t)xl     ^  Q 

uniformly  in  t  as  n  -»•  «. 

From  lemma  5.4   (ii)   and  (5.67),   it  is  clear  that  the  right 
hand  side  of  (5.65)   converges  to  0  as  n  ^  »  uniformly  in  t. 
Thus,   from  (5.65),  we  conclude  that 

|p^(t)x-P( t)x|^  ^  0  as  n  >  »  uniformly  in  t  e  [0,T]. 
That  is,   P^(t)  ^  P(t)   strongly  in  l^(Q)   as  n  >  »  uniformly 
in  t.     From  the  Riccati  integral  equation  in  theorem  5.1,  it 
can  be  easily  seen  that  for  x,  y  e  L^(Q), 
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(P(t)x,y)^  =  (S(s-t)x,FS(s-t)y)^ds 

-  lim         (BP   (s)S(s-t)x,BP^(s)S(s-t)y)  ds 
+  (S(T-t)x,GS(T-t)y )q       uniformly  in  t. 

This  completes  the  proof  of  theorem  5.3. 
Finally,  we  are  in  a  position  to  complete  the  proof  of 
theorem  4.1  without  (A.6). 

Proof  of  theorem  4.1   (ii).     From  theorem  4.1   (ii)  under 

2 

(A.6),  we  obtain  that  for  x,   y  €  L  (Q),   t  €  [0,T], 

(5.68)  (P^(t)x,y)2  =  !l^^^^^rt)x,F^$^(s,t)y)^65 

+  jl  (BP   (s)$  (s,t)x,BP^(s)$^(s,t)y)„ds 
'  t        n        n  n        n  i 

+  ($^(T,t)x,G^$^(T,t)y)^. 
From  theorem  5.3,  we  have 

(5.69)  (P^(t)x,y)Q  ^   (P(t)x,y)Q  as  n  ^  »  uniformly  in  t. 

By  lemma  5.4   (ii),   it  is  clear  that 

(5.70)  ($^(s,t)x,F^$j^(s,t)y)2ds 

/'J  ($(s,t)x,F$(s,t)y)2ds  as  n  CO, 

and   ($^'T,t)x,G^$^(T,t)y)2  >   ( $ ( T , t ) x ,G$ ( T , t ) y ) ^  as  n  >  «. 
From  (4.17),  we  recall  that  for  x  €  L^(Q),   t  €  [0,T), 
u°(»,t;x)   =  BP(»)$(»,t)x  and  u°(«,t;x)    =  BP^  (  •  )      (  •  ,  t )  x  . 
By  lemma  5.4   (i),   it  is  easily  seen  that 
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(5.71)  (u°(.,t;x),u°(.,t;Y))j^^^j^^ 

-»•   (u°(  •  ,  t;x)  ,u°(  •  ,  t;Y)  )  j^^^j       as  n       =  . 

In  fact, 

|(u°(.,t;x),u°(.,t;y))^^^^^^^-(u°(.,t;x),u°{.,t;y))  ^^^^j^J 

<  |(Un^-^t;x),u°(.,t;y))^^^^j^^-(u°(.,t;x),u°(.,t;y))^^^^j^^| 

+    |(u°(.,t;x),u°(.,t;y))j^^^^^^-   (  u°(  •  ,  t ;  x  )  ,  u°(  •  ,  t ;  y  )  )  ^  ^  ^ 

=  |(u°(.,t;x),u°(.,t;y)-u°(.,t;y))^^^^j^^| 

+  |(u°(.,t;x)-u°(.,t;x),u°(.,t;y))^^^^j^^| 

->■  0  as  n  ->•  "  which  follows  from  lemma  5.4  (i). 

Thus,  we  may  rewrite   (5.71)  as 

(5.72)  (BP^(s)$^(s,t)x,BP^(s)$^(s,t)y)^ds 

■»•         (BP(s)$(s,t)x,BP(s)$(s,t)y)   ds  as  n  «. 
Hence,  by  taking  limits  in  (5.68),  we  obtain 

(P(t)x,y)^   =  ($(s,t)x,F$(s,t)y)^ds 
T 

+         (BP(s)$(s,t)x,BP(s)$(s,t)y)^ds  +  ($(T,t)x,G$(T,t)y)^. 

This  completes  the  proof  of  theorem  4.1  (ii). 

From  theorem  4.1   (ii),  we  obtain  the  following 
corollary  without  any  smoothness  assumptions  on  F  and  G. 
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Corollary .     For  each  t  €   [0,T),  we  have 

(i)  P(t)    is  self-adjoint  and  positive-definite, 

(ii)  the  minimal  value  of  the  functional  cost  J(.  on  the 
interval   [t,T]   with  initial  data  x  is  given  by 
J^(u°(. ,t;x) ,y°( . ,t;x) )   =^  (P(t)x,x)g. 

VJe  give  examples  of  how  to  obtain  sequences  of 

smoothing  operators  F^j  and  G^^  converging  to  F  and  G  strongly 
2 

in  L  (a)  respectively. 

Example  1.     Since  the  operators  A  and  A*  generate 

2 

strongly  continuous  semigroups  on  L  (Q),   there  exists  a 
number         such  that  {r   |    r  >  u)^}    is  contained  in  the 

resolvent  sets  of  A  and  A*   [1].     Thus,   the  operator 

—  1  —  1  2 

(n+(jj^-A*)      (n+uj^-A)        is  bounded  on  L   (Q)    for  each  positive 

integer  n. 

Let         =   (n+u)^)^(n+a)^-A*)~-'-(n+a)Q-A)~-'-   for  each  n.     Then  it 

is  clear  that  R^^  is  self-adjoint  and  positive-definite  for 

2 

each  n.     Moreover,   R^  maps  L   (Q)    into  D(A*)   for  each  n. 

On  the  other  hand,   it  is  well  known   [1]  that 

2 

R    >  I  strongly  in  L  (Q)   as  n  ^  « .   Let  F     =  R  F  and  G     =  R  G 

n         n  n  n 

for  each  n.     Then  F^^  and  G^^  are  self-adjoint,  positive- 
definite  and  bounded  operators  on  L^(Q)  which  map  L^{Q)  into 
D(A*)   for  each  n. 
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2 

For  each  x  €  L  (Q),  we  have 


F^x-Fx|^  =    |R^Fx-Fx|^       0  and 
Ig^x-GxL  =   |r  Gx-GxL  ■>  0  as  n  •>  00. 

This  implies  that  { F^}  ^  and  {G^^}^  satisfy  (5.53). 

Example  2.     Let        =  A  with  zero  boundary  condition  on 

m    ^  2  2  1 

Q,  where  A  =     Z  ^— =■    with  domain  D(A,)   =  H   (Q)  H  H  (Q). 

j=l  5x^  ^  ° 

Then  it  is  clear  that  -A^  is  self-adjoint  and  positive- 

1 
2 

definite.     Thus,  -A,     is  also  self-adjoint  and  positive- 

i 

2  1 

definite  with  domain  D(-A^)  =  H^(Q). 

1 
2 

It  is  well  known  that  A^  generates  a  strongly 

2 

continuous  semigroup  on  L  (Q).     Thus,   for  some  real  number 

0)^,  the  set  (r  I   r  >  u^}   is  contained  in  the  resolvent  set 

1  1 
of  A^.     That  is,   ("+Wo"'^d^~^        bounded  on  L^(Q)   for  each 

1 

positive  integer  n.     Let        =  ( n+w^)  (  n+w^^-A^ ) for  integer 

I 
2 

n  >  1,     Then,  since  -A^  is  self-adjoint  and  positive- 
definite,  so  is        for  each  n  >  0. 

2  1 
Moreover,   R^  maps  L  (Q )    into  H^(Q)   for  each  n. 

It  is  well  known  [1]    that  I  strongly  in  L^(Q)  as 

n      CO.     Let  F^  =  R^F  and  G^  =  R^G  for  each  n.     As  in  example 

1,   it  is  easily  seen  that  {f^}^  and  {G^}^  satisfy  (5.53). 
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